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Executive Summary 



We propose an experiment (Mu3e) to search for 
the lepton flavour violating (LFV) decay /i + —¥ 
e + e~e + . We aim for an ultimate sensitivity of one 
in 10 16 /i-decays, four orders of magnitude better 
than previous searches. This sensitivity is made 
possible by exploiting modern silicon pixel detect- 
ors providing high spatial resolution and hodo- 
scopes using scintillating fibres and tiles providing 
precise timing information at high particle rates. 

Existing beamlines available at PSI providing 
rates of order 10 8 muons per second allow to test 
for the decay /i + —¥ e + e~e + in one of 10 15 muon 
decays. In a first phase of the experiment, we 
plan to make use of this and establish the experi- 
mental technique whilst at the same time pushing 
the sensitivity by three orders of magnitude. 

The installation of a new muon bcamlinc at the 
spallation neutron source is currently under dis- 
cussion at PSI. Such a High Intensity Muon Beam 
(HiMB) will provide intensities in excess of 10 9 
muons per second, which in turn are required to 
reach the aimed sensitivity of B(/i + — > e + e~e + ) ~ 
lCr 16 . 

The proposed experiment is highly complement- 
ary to other LFV searches for physics beyond the 
standard model, i.e. direct searches performed at 
the Large Hadron Collider (LHC) and indirected 
searches in the decay of taus and muons, such 
as the decay /i + — > e + 7, which is the subject 
of the MEG experiment currently in operation 
at PSI. The proposed experiment for the search 
/i + — > e + e~e + will test lepton flavour violating 
models of physics beyond the Standard Model 
with unprecedented sensitivity. 

This sensitivity is experimentally achieved by a 
novel experimental design exploiting silicon pixel 
detectors based on High Voltage Monolithic Act- 



ive Pixel Sensors (HV-MAPS). This technology 
provides high granularity, important for precision 

tracking and vertexing, and allows one to signi- 
ficantly reduce the material budget by thinning 
down the sensors and by integrating the hit digit- 
isation and readout circuitry in the sensor itself. 
The detector geometry is optimized to reach the 
highest possible momentum resolution in a mul- 
tiple Coulomb scattering environment, which is 
needed to suppress the dominating background 
from the radiative muon decay with internal con- 
version, fi — > eeevv. The time information of the 
decay electron^] obtained from the pixel detector 
is further improved by a time-of-flight system con- 
sisting of a scintillating fiber hodoscope and tiles 
with Silicon Photo-Multipliers (SiPM) for light 
detection. By combining both detector systems 
accidental background can be reduced below the 
aimed sensitivity of B(/i + — > e + e~e + ) ~ 1CP 16 . 

We will complete the sensor development and 
start constructing the detector in 2013, in order 
to be ready for first exploratory data taking at an 
existing beam line with a first minimal detector 
setup in 2015. A detector capable of taking data 
rates of order 10 8 muons per second and capable 
of reaching a sensitivity of B(/i + — > e + e~e + ) ~ 
10~ 15 will be available in 2016. This Phase I de- 
tector is the main focus of this proposal. 

In Phase II, beyond 2017, the experiment will 
reach the ultimate sensitivity by exploiting a pos- 
sible new high intensity muon beamline with an 
intensity of > 2 ■ 10 9 muons per second. In the 
absence of a signal, LFV muon decays can then 
be excluded for B(>+ e+ e -e+)< 10~ 16 at 90% 
confidence level. 

1 Here and in the following, the term "electron" denotes 
generically both decay electrons and positrons. 
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Part I 



Chapter 1 




Motivation 



In the Standard Model (SM) of elementary 
particle physics, the number of leptons of each 
family (lepton flavour) is conserved at tree level. 
In the neutrino sector, lepton flavour violation 
(LFV) has however been observed in the form 
of neutrino mixing by the Super-Kamiokande \U , 
SNO |], KamLAND [3] and subsequent experi- 
ments. Consequently, lepton flavour is a broken 
symmetry, the standard model has to be adap- 
ted to incorporate massive neutrinos and lepton 
flavour violation is also expected in the charged 
lepton sector. The exact mechanism and size of 
LFV being unknown, its study is of large interest, 
as it is linked to neutrino mass generation, CP 
violation and new physics beyond the SM (BSM) . 

The non-observation of LFV of charged leptons 
in past and present experiments might at a first 
glance be surprising, as the mixing angles in the 
neutrino matrix have been measured to be large 
(maximal) . This huge suppression of LFV effects 
is however accidental and due to the fact that (a) 
neutrinos are so much lighter than charged leptons 
and (b) the mass differences between neutrinos 
(more precisely the square of the mass differences) 
are very small compared to the W-boson mass. 

The situation completely changes if new 
particles beyond the SM are introduced. If e.g. 
SUSY is realized at the electroweak scale, the 
scalar partners of the charged leptons (sleptons) 
will have large masses, and if not fully degener- 
ate, induce LFV interactions via loop corrections. 
These LFV effects from new particles at the TeV 
scale arc naturally generated in many models and 
are therefore considered to be a prominent signa- 
ture for new physics. 



In many extensions of the SM, such as grand 
unified models |4}{6], supersymmetric models 7 



(see Figure 2.2 1, left-right symmetric models [8 
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[10], models with an extended Higgs sector 
and models where electroweak symmetry is broken 
dynamically 12 , an experimentally accessible 
amount of LFV is predicted in a large region of 
the parameter space. 

Seesaw and Left-Right symmetric (supersym- 
metric) models are good candidates for realising 
grand unification, which also unify quark and 
lepton mass matrices. Moreover, it has been 
shown that LFV effects in the low energy limit 
can be related to mixing parameters at the GUT 
scale or to heavy Majorana masses in these mod- 
els 



17 18 



Seesaw models are therefore very at- 
tractive in the context of LFV as they are also able 
to naturally explain the smallness of the masses of 
the left handed neutrinos. In this context the re- 
cent results from neutrino oscillation experiments 
are very interesting, as they measured a large mix- 
ing angle 6*13, which enhances the LFV-muon de- 
cays in most models which try to explain the small 
neutrino masses and the large mixing. 

Currently the most accurate measurement is 
provided by t he Daya Bay reactor neutrino ex- 
periment [19] yielding sin 2 (26» 13 ) = 0.089 ± 
O.OlO(stat) ± 0.005(syst), excluding the no- 
oscillation hypothesis at 7.7 standard deviantions. 
The Daya Bay measurement is in good agreement 
with measurements by the RENO [20], Double 
Chooz [21] and T2K [22] experiments. These 



results are very encouraging, as large values of 
sin 2 (26*i 3 ) lead to large LFV effects in many BSM 
models. 
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Decay channel Experiment Branching ratio limit Reference 

fi^e-/ MEGA <1.2-Kr n [l3 

MEG <2.4-10~ 12 14 

/i^eee SINDRUM < 1.0 • 10~ 12 TE 

fiAu-^eAu SINDRUM II < 7 ■ 1CT 13 T5 



Table 1.1: Experimental limits on LFV muon decays 



The observation of LFV in the charged lepton 
sector would be a sign for new physics, possibly 
at scales far beyond the reach of direct observa- 
tion at the large hadron collider (LHC). Several 
experiments have been performed or are in opera- 
tion to search for LFV in the decays of muons or 
taus. Most prominent are the search for the radi- 
ative muon decay fi — > e~f 13p4{|23p4 , the decay 



grade the experiment to further improve the sens- 
itivity are currently under discussion. 

In the near future the DeeMe experiment 
[44] at J-PARC plans to improve the current 
muon-to-electron conversion limit of B(/x Au 
eAu)<7-lQ- 



-13 



1G 



15 , the conversion of captured muons to by COMET at J-PARC 45 



fi — S- eee 

electrons 16 and LFV tau decays |25}|43 



by almost two orders of mag- 
nitude. By the end of the decade this limit could 
be improved by even four orders of magnitude 

and Mu2e at Fer- 



milab 46 47 



The recent search performed by the MEG 

-12 14 



Collaboration yields B(n ->• 67) < 2.4 • 10 
and sets currently the most stringent limit on 
many LFV models. The MEG collaboration plans 
to continue operation into 2013 in order to in- 
crease the number of stopped muons and to reach 
a sensitivity of a few times 10~ 13 . Plans to up- 



Selected limits for lepton flavour violating muon 
decays and muon-to-electron conversion experi- 
ments, which are of high relevance for the pro- 
posed experiment, are shown in Table |1.1| A 
search for the LFV decay \i — > eee with an unpre- 
cedented sensitivity of < 10~ 16 as proposed here 
would provide a unique opportunity for discover- 
ies of physics beyond the SM in the coming years. 
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Chapter 2 




Theory 



In the SM, charged lepton flavour violating reac- 
tions are forbidden at tree level and can only be in- 
duced by lepton mixing through higher order loop 
diagrams. However, the dominant neutrino mix- 
ing loop diagram, see Figure |2.1[ is strongly sup- 
pressed in the SM with B <§; 10~ 50 and thus giv- 
ing potentially high sensitivity to LFV reactions 
in models beyond the Standard Model (BSM). 

Such an example is shown in Figure [2T2] where 
a 7/Z-penguin diagram is shown with new su- 
persymmetric (SUSY) particles running in a loop. 
These loop contributions are important basically 
for all models, where new particle couplings to 
electrons and muons are introduced. Lepton flavor 
violation can also be mediated by tree couplings 
as shown in Figure |2.3[ These couplings could 



be mediated by new particles, like Higgs particles 
or doubly charged Higgs particles, R-parity viol- 
ating scalar neutrinos or new heavy vector bo- 
sons, the latter being motivated by models with 
extra dimensions |48j|49] . These models usually 



yVZ, 



+ 



X 



Figure 2.2: Diagram for lepton flavour violation 
involving supersymmetric particles. 



also predict semihadronic decays of tau leptons or 
the muon conversion process fiq —> eq', which is 
experimentally best tested in muon capture exper- 
iments. 




X 



Figure 2.1: Feynman diagram for the fi — > eee pro- 
cess via neutrino mixing (indicated by the cross). 




Figure 2.3: Diagram for lepton flavour violation 
at tree level. 
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The lepton flavor violating three electron decay 
of the muon can be mediated, depending on the 



model, via virtual loop (Figure 2.2) and box dia- 
grams or via tree diagrams (Figure 2.3 1. The most 
general Lagrangian for this decay can be paramet- 
erized as [50 Q 



2 L 

+ gi (JTrcl) (e]ie L ) 
+ 92 (~pLe R ) (eZe R ) 
+ 93 (7AR7 M e fl ) (e^e^) 
+ 34 (Ml7 Me i) (ei7M e i) 
+ 35 (Wr^cr) (eLl^e L ) 
+ 96 {WI^sl) (eR-j^eR) + H.c] 

(2.1) 

The form factors Ar : l describe tensor type (di- 
pole) couplings, mostly acquiring contributions 
from the photon penguin diagram, whereas the 
scalar-type (171,2) and vector-type (173 — ge ) form 
factors can be regarded as four fermion contact 
interactions, to which the tree diagram contrib- 
utes in leading order. In addition also off shell 
form factors from the penguin diagrams, which 
are not testable in the fj, — > ej decay contribute 
In case of non-zero dipole and four- 



51 



to 31-36 

fermion couplings also interference effects have to 
be considered, which can be exploited to investig- 
ate violation of time reversal (T-invariance) . 

By neglecting higher order terms in m e , the 
total branching ratio of the decay can be expressed 
by: 



B(fi — > eee) 



9l 



9l 



m,, 

32 eA 2 (In -£ 

mi 



+ 2 (33 + 34) + 35 + 3e 
11/4) 



+16 t] eA 
+8 rf eA 



9l 



9l 



9l 



9l 



(2.2) 



where the definition A 2 = A\ 



A\ is used. The 



term proportional to A 2 is logarithmically en- 
hanced and can be assigned to the photon penguin 



1 A representation of Lagrangian containing explicitly 
the contributions from the loop and box diagrams can be 
found in |51|. 



diagram. The constants 77 and 77' are T-violating 
mixing parameters. In case of a signal, the dif- 
ferent terms can be measured from the angular 
distribution of /1 — > eee decay particles using a 
polarized muon beam. 



2.1 Comparison fi 



eee versus 



In the decay [i — > ey physics beyond the SM is 
only tested by photon penguin diagrams, in con- 
trast to fi — > eee where also tree, Z-penguin and 
box diagrams contribute. To compare the new 
physics mass scale reach between the processes 
/i — > eee and fi — > e^j a simplified model is chosen; 
it is assumed that the photon penguin diagram 
Figure |2.2| and the tree diagram Figure |2.3| are 
the only relevant contributions. The Lagrangian 
then simplifies tcP] 



J LFV 



^r^^lF^ 



7— pcngum 



(k + 1)A 2 



{UlI^cl) (e^e^) 



tree 

(2.3) 



where for the contact interaction ("tree") term ex- 
cmplarily a left-left vector coupling is chosen. In 
this definition a common mass scale A is intro- 
duced and the parameter k describes the ratio of 
the amplitudes of the vector-type (tree) term over 
the tensor (7 — penguin) term. 

Limits on the common mass scale A as ob- 
tained from the experimental bounds on B(/i — > 
e 7 )< 2.4- 1CT 12 (90% CL MEG 2011) andBQu -> 
eee)< 1.0 • 10" 12 (90% CL SINDRUM) arc shown 
in Figure [2~4] as function of the parameter k. Ex- 
perimentally, for small values of n (dipole coup- 
ling) the mass scale A is best constrained by the 
MEG experiment whereas the four fermion con- 
tact interaction region with n > 10 is best con- 
strained by the SINDRUM experiment. 

For comparison also a hypothetical ten times 
improved limit is shown for the MEG experiment 
(post-upgrade) and compared to the sensitivities 
of the proposed [i — > eee experiment of 10~ 15 
(phase I) and 10~ 16 (phase II). It can be seen 
that in this simple model comparison high mass 
scales A will be best constrained by the proposed 



J A similar study was presented in 52 
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10 10 
10* 





Figure 2.4: Experimental limits and projected lim- 
its on the LFV mass scale A as a function of the 



Figure 2.5: Experimental limits and projected lim- 
its on the LFV mass scale A as a function of the 



parameter re (see equation 2.3 1 assuming negligible 
contributions from Z° penguins; based on [52] . 

/i — > eee experiment for all values of re already in 
phase I. 

In case of dominating tensor couplings (A ^ 
0, re — > 0) a quasi model independent relation 
between the /i — > eee decay rate and the /i — > e^ 
decay rate can be derived: 



B(/i — > eee) 
B(fi -> erf) 



0.006 



(2.4) 



This ratio applies for many supersymmetric mod- 
els, where LFV effects are predominantly medi- 
ated by gauge bosons and where the masses of 
the scalar leptons or gauginos are of electroweak 
scale. In these models, which are already heavily 
constrained or even excluded by the recent LHC 
results, the sensitivity of the proposed MuSe ex- 
periment in terms of branching ratio has to be two 
orders of magnitude higher than that of the MEG 
experiment in order to be competitive. 

2.1.1 Z-penguin Contribution 

However, besides the tree and 7-penguin dia- 
grams also the Z-penguin diagram can signific- 
antly contribute to the process fi —> eee. The 
Z-penguin diagram is of particular importance if 
the new physics scale is higher than the electro- 
magnetic scale, as can be easily derived from a 
dimensional analysis. The enhancement of the 
Z-pcnguin contribution over the 7-penguin con- 
tribution and its non-decoupling behaviour when 



parameter re (see equation 2.3 1 assuming contri- 
butions from Z° penguins ten times larger than 
the photon contribution. 



going to high mass scales was discussed for Little 



Higgs models 53 54 as well as for several SUSY 
models 55-59). SUSY models with i?-parity viol- 
ation and right handed neutrinos received recently 
quite some attention in this context, as approxim- 
ate cancellations of different Z-diagram contribu- 
tions are not present in extended Minimum Super- 
Symmetric Standard Models (MSSM). 

The effect of such an enhanced Z-penguin coup- 
ling, where the LFV contribution to the fj, — > eee 
process is exemplarily enhanced by a factor of 
ten relative to the photon-penguin contribution, 
is shown in Figure |2.5| It can be seen that the 
sensitivity of the /1 — > eee process to new physics 
is significantly enhanced at small values of re and 
that in such sensitivity of 10 14 is already 

sufficient to be competitive with the fj, ^ ej pro- 
cess with a sensitivity of a few times 10~ 13 . 

2.2 Discussion of Specific Models 

In the following, selected models are discussed in 
more detail in the context of the proposed exper- 
iment. 

2.2.1 Inverse Seesaw SUSY Model 

Despite the fact that the most simple supersym- 
metric models with light squarks and gluinos were 
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1500 2000 2500 3000 



m o = Mi/2 [GeV] 

Figure 2.6: Inverse Seesaw SUSY Model: Contri- 
butions to -B(/i — > eee) as a function of mo = M1/2 
for a degenerate singlet spectrum with Mr = 
lOTeV and M = 10 11 GeV. The rest of the 
cMSSM parameters are set to A — —300 GeV, 
B = 0, tan/3 = 10 and sign(^i) = +. Solid lines 
represent individual contributions, 7 (black), Z 
(blue) and h (red) whereas the dashed lines rep- 
resent interference terms, 7 — Z (green), 7 — h 
(purple) and Z — h (orange). Note that in this 
case h includes both Higgs and box contributions. 



Taken from 58 



60 84 



recently excluded by LHC experiments 
supersymmetry can still exist in nature, just at 
higher mass scales or in more complex realisations. 
In many of these realisations with a non-minimal 
particle content the Z-penguin contribution dis- 
cussed above gets significantly enhanced. 

As a first example results obtained by a super- 
symmetric model with an inverse seesaw mechan- 
are discussed here 



ism 



model 



58 



85 



The inverse seesaw 
constitutes a very appealing alternat- 



ive to the standard seesaw realization and can be 
embedded in a minimal extension of the MSSM 
by the addition of two extra gauge singlet su- 
perfields, with opposite lepton numbers. Similar 
to other models, e.g. flavour violating Higgs de- 
cays in the MSSM, the Z-penguin exhibits here 
a non-decoupling behaviour, which is shown in 
Figure |2.6| for an effective right-handed neut- 
rino mass of M — 10 11 GeV and degenerate 
sterile neutrino masses of Mr = 10 TeV. At 
small mass scales mo = M1/2 of the constrained 
MSSM (cMSSM) the photon-mediated penguin 
contribution clearly dominates over the other con- 
tributions from Higgs-mediated penguin and Z- 
mediated penguin diagrams. This picture com- 
pletely changes at higher mass scales above 200- 



T 
a. 

« 10 



10- 9 

10 -io 

10-" 
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500 1000 1500 2000 2500 3000 3500 4000 



Figure 2.7: Supersymmetric SU(3) C x SU(2) L x 
U(1)b-l x U(1)ji Model: Branching ratios of 
lepton flavour violating processes as a function 
of m vR for m = 800 GeV, Af 1/2 = 1200 GeV, 
tan/3 = 10, A = 0, v R =10 TeV, tan/3 fl = 1.05, 
fi R = -500 GeV, m Ali = 1000 GeV. The dashed 
red line is the predicted branching ratio for fi — > 
e7 and the dashed blue line for /1 — > eee. Taken 
from |59| . 



300 GeV, where the Z-mediated penguin diagram 
becomes dominant. The non-decoupling beha- 
viour of the Z-penguin is clearly visible which will 
allow to test this model at any SUSY mass scale 
for the seesaw parameters given in this example 
at phase II of the proposed experiment. 

2.2.2 Supersymmetric SU(3) c x SU(2) l x 
U(1) B -l x U(1) R Model 

This model represents a supersymmetric version 
of the SM, minimally extended by additional 



U(1)b-l x U(1)r symmetry groups 86 87 . This 



model includes the generation of light neutrino 
masses by the seesaw mechanism, can explain the 
observed large neutrino mixing angles and can 
be easily embedded into a 5*0(10) grand unified 
theory. This model predicts an additional light 
Higgs particle, which is expected to mix with the 
lightest MSSM Higgs particle, and has been re- 
cently studied also in the context of LFV pro- 
Also in this study it is found that at 



59 



cesses 

high SUSY mass scales the photon-mediated LFV 
penguin diagrams are more suppressed than the 
Z-mediated LFV penguin diagrams and that this 
suppression scales with m| / m susY as naively 
expected from a dimensional analysis. Branch- 
ing ratio predictions for the processes /x — > eee 
and /i — >■ e7 are shown in Figure |2.7| as func- 
tion of the right-handed neutrino mass m UR for 
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the SUSY model parameters as given in the figure 
caption. Also here it can be seen that for high 
masses m VR > 300 GeV the Z-mediated penguin 
diagram starts to contribute dominantly to the 
/i — > eee process and that for m UR > 1000 GeV 
the \i — > eee process is expected to have an even 
higher branching fraction than \i — > e'y. For even 
higher masses the non-decoupling behaviour is vis- 
ible in the /i — > eee prediction. 

2.2.3 Other Models 

The above discussed enhancement of the Z- 
mediated penguin diagram appears also in Little 
Higgs Models with T-parity (LHT) where ratios 
B(/x — > eee)/B(// —> e^y) « 0.02 — 1 have been pre- 
dicted [53j[54)[88)[89], or in Left- Right Symmetric 
models with additional Higgs triplets. LFV inter- 
actions in Higgs-triplet models can be also gener- 
ated directly in tree diagrams, see Figure [273] 

these LFV violating effects are studied 
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in a model where the Higgs triplet is responsible 
for neutrino mass generation. Figure |2.8| shows 
the predicted branching ratios for each of the three 
LFV muon processes and for different realisations 
of the neutrino mass hierarchy. Note that the ab- 
solute value of the branching ratios depends on 
the mass scale M and can vary. For the hier- 



archical case, Figure 2.8 b), all branching ratios 



are expected to be of similar size whereas for the 



degenerate, Figure 2.8 a), and the inverted case, 
Figure 2.8 c), the fi — > eee branching ratio dom- 
inates in the allowed region of U e s. As the LFV- 
mediating Higgs triplet boson does not couple to 
quarks, the /x — > eee decay is enhanced compared 
to the ji — > ej decay and the muon-to-electron 
conversion processes, which are both loop sup- 
pressed. 



This enhancement of the LFV tree diagram is 
also found to be large in extra dimension mod- 
els 48p0 or models with new heavy Z bosons. In 
Randall-Sundrum (RS) models |48 , flavor chan- 
ging neutral currents (FCNCs) arise already at the 
tree level. This is caused by the flavor-dependent 
couplings of these gauge bosons, due to their non- 
trivial profiles in the extra dimension. Moreover, 
FCNCs arise through the exchange of the Higgs 
boson, as due to the contribution to the fermion 
masses from compactification, there is a misalign- 
ment between the masses and the Yukawa coup- 
lings. 

Electroweak precision observables suggest that 
for RS models featuring the Standard Model 
gauge group, the new-physics mass scale Mkk 
(the scale of the Kaluza Klein excitations) should 



not be lower than O(10 TeV) at 99% CL El -93 



Thus, without additional structure/symmetries, 
the experimental situation suggests that it could 
be challenging to find direct signals from RS mod- 
els at the LHC. In such a situation, precision 
experiments, like the measurement of the decay 
/i — > eee, will furnish the only possibility to see 
the impact of warped extra dimensions. 

2.3 Theory Summary 

The search for the decay /i — > eee is in itself of fun- 
damental interest and might reveal surprises not 
foreseen by any of the models discussed above. 
This search is largely complementary to other 
LFV searches, in particular to the decay \i — > ej 
and to the /x — > e conversion in muon capture 
experiments. In a wide range of models for phys- 
ics beyond the standard model, highest sensitiv- 
ity in terms of branching ratio is expected for the 
/i — > eee decay process. 
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Figure 2.8: Branching ratios of B(/i — > eee), B(// — > ej) and muon conversion B(/i — > e) in different Higgs 
triplet scenarios with a degenerated, hierarchical or inverted mass hierarchy of the neutrinos as function 
of the neutrino mixing matrix element U e 3 and for the model parameters: M = 200 GeV, A = 25 eV, 
and m v = 0.1 eV for the degenerate case. These plots were taken from [TT] and the U e ^ constraints 



(green bands) as obtained from 19 were added posterior 
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Chapter 3 




Experimental Situation 



3.1 SIN DRUM Experiment 

The SIN DRUM experiment was in operation at 
PSI from 1983-86 to search for the process \i — > 
eee. No signal was found and the limit B(/i — > 
eee)< 10" 12 was set at 90% CL [15], assuming a 
decay model with a constant matrix element. 

The main components of the experiment were 
a hollow double-cone shaped target of dimension 
58 x 220 mm 2 to stop surface muons of 28MeV/c 
in a solcnoidal magnetic field of 0.33 T, five layers 
of multiwire proportional chambers and a trigger 
hodoscope. The main tracking parameters which 
were most relevant for the search sensitivity of the 
experiment are shown in Table [3~T] 

The time resolution obtained by the hodoscope 
of less than 1 ns was, together with the achieved 
momentum resolution, sufficient to suppress the 
accidental background completely. 

After all selection cuts, no candidate event was 
seen by the SINDRUM experiment. The sensitiv- 
ity of the experiment was mainly determined by 
the \i — > eeevv background process and estimated 
as 5 • 10 -14 94 ; the obtained limit was basically 
given by the limited number of muon stops. 



In the first running period in the year 2008 



3.2 MEG Experiment 

The MEG experiment at PSI is in operation since 
2008 and is searching for the LFV decay /i — > ej. 
The main components used for event reconstruc- 
tion are drift chambers for positron detection and 
a liquid xenon calorimeter for photon detection. 
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about 10 14 muons were stopped on target 
No signal was found and a limit on the decay of 
B(/x -> e 7 )< 2.8 • 10" 11 (90 % CL.) was set. 

After upgrading the detector the search sensit- 
ivity and the limit was improved using data taken 
in the years 2009/2010 to B(> -> e 7 )< 2.4 • 10" 12 
(90% CL.) [14]. 

The dominant background contribution for [i — > 
e 7 comes from accidentals where a high energy 
photon from a radiative muon decay or from 
a bremsstrahlung process is recorded, overlayed 
with a positron from the upper edge of the Michel 
spectrum. This accidental background mainly de- 
termines the final sensitivity of the experiment. 

The amount of background is predominantly de- 
termined by the timing, tracking and energy resol- 
ution. Selected resolution parameters as achieved 
in the 2009 run are summarized in Table 13.21 
The MEG experiment will continue operation un- 
til middle of 2013. The final sensitivity is expec- 
ted to be a few times 10~ 13 . The collaboration 
has started to discuss possible upgrades to fur- 
ther improve the sensitivity by about one order of 
magnitude. These numbers are to be compared to 
the bound from the earlier MEGA experiment of 
B(/x -> e 7 )< 1.2 



10" 



13 



The study of the /i — > e 7 decay sets stringent 
bounds on models predicting new heavy particles 
mediating LFV dipole couplings. These dipole 
couplings can also be tested in the process /z — s- 
eee, where the sensitivity is reduced by a factor 
of ^ (In (m 2 /m 2 ) - 11/4) = 0.006 (note however 
that for /i — > eee also box diagrams, Z°-mediated 
penguin diagrams and tree digrams contribute as 
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SINDRUM PARAMETER 



Value 



rcl. momentum resolution <J p /p 5.1% (p = 50MeV/c) 

rcl. momentum resolution a p jp 3.6% (p = 20MeV/c) 

polar angle <jg 28mrad (p = 20MeV/c) 

vertex resolution <Jdca ~ 1 mm 

MWPC layer radiation length in X 0.08 % - 0.17 % 



Table 3.1: SINDRUM tracking parameters taken from 15 . 



described in chapter [2| . In the case that the LFV 
dipole couplings are dominant, the projected sens- 
itivity of 10 -13 of the MEG experiment corres- 
ponds accordingly to a sensitivity of about 10~ 15 
in the search for the [i — > eee decay and the en- 
visaged sensitivity of B(fj, — >• eee) = 10 -16 corres- 
ponds to more than one order of magnitude higher 
sensitivity compared to the MEG experiment. 



3.3 Muon Conversion Experiments 

Muon to electron conversion experiments /i — > e 
on nuclei exploit the clear signature of monochro- 
matic electrons. Differently to the search for LFV 
muon decays, which are performed using DC anti- 
union beams in order to reduce accidental back- 
grounds, muon conversion experiments are per- 
formed using pulsed muon beams to reduce the 
rapidly decaying pion background. A limitation 
of this type of experiment is the background from 
ordinary decays of captured muons with large nuc- 
lear recoil and from pions. 

The most stringent limits for muon-electron 
conversion on various nuclei have been obtained 
by the SINDRUM II collaboration |B|97|98] . The 



MEG PARAMETER 2011 PUBL. 


Value 


rel. momentum resolution a p /p 


0.7% (core) 


polar angle ag 


9mrad 


azimuthal angle 


7mrad 


radial vertex resolution ctr 


1.1 mm 


long, vertex resolution : <iz 


1.5 mm 



Table 3.2: Best MEG tracking parameter resol- 
utions achieved in the year 2009/2010. The res- 
olutions are given for positrons of 53MeV/c mo- 
mentum. Values taken from 



strongest limit has been set using a gold target 
B(/i Au^eAu)<7- 10" 13 [l6 . 

Similar to the /i — > eee process, the sensitivity 
to dipole couplings in muon conversion is reduced 
by about a em compared to the more direct /i — > 
ej search. However, new experiments planned 
at Fermilab (Mu 2e |i6| [99lp0]) and at J-PARC 
{COMET |45|101||102| and PRISM [IMpuI] ) aim 
for branching ratios of 10 -16 or smaller relative 
to the captured muon decay and have a higher 
sensitivity to LFV dipole couplings than the run- 
ning MEG experiment. Similar to the fi — > eee 
process, also four-fermion couplings are tested in 
fj, —> e conversion experiments. These couplings 
involve light quarks and are thus complementary 
to all other LFV search experiments. 

The Mu2e and COMET experiments are ambi- 
tious projects and are expected to come into oper- 
ation at earliest by the end of this decade. In a few 
years time the DeeMe experiment at J-PARC 44 
will start taking data, aiming for a sensitivity for 
muon-to-electron conversions of 10 -14 . 



At Osaka university, the MuSIC project 105 
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106 aims for a very high intensity DC muon beam 
using a high-field capture solenoid around a thick 
conversion target. One of many possible users of 
that beam is a \x — > eee experiment. However, an 
experimental concept has not yet been presented. 



3.4 LFV in r Decays 

A wide variety of LFV decay channels are open in 
r decays. These decay modes have been extens- 
ively explored at the B-factories, producing limits 
on branching ratios of a few 10~ 8 , see Table 3.3 
and Figure |3.1| The next generation of B exper- 
iments at e + e~ colliders could push these limits 
down by one to two orders of magnitude. For cer- 
tain channels such as r — > fj-fifj,, the LHCb exper- 
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Figure 3.1: Limits on LFV r decays. Taken from 95 



imcnt could also be competitive given the lumin- 



osity expected in the coming years 107 108 



3.5 LFV at the Large Hadron Col- 
lider 

LFV signatures might be observed at the LHC 
if e.g. supersymmctric particles are discovered, 
which naturally generate LFV couplings in slepton 
mass mixing. Consequently, if sleptons arc light 
enough to be produced in pairs, different lepton 
flavors might show up in decay chains such as: 
£+t -> £+£-'x a X°- 

Known and new scalar or vector particles could 
also have lepton violating tree couplings and might 
be directly reconstructed from resonance peaks: 
H -> 11' or Z' -4 11'. Due to the existing 
bounds on flavor changing processes, these LFV 
decays are small and difficult to detect above 



the large background from VFVF-production with 
subsequent leptonic decays. It seems however, 
that with high enough luminosities, the LHC can 



e.g. go be yond the LEP bounds |109^112] on LFV 
Z decays |113| . 

If new particles exist at the TeV mass scale, 
i.e. in the discovery reach of the LHC, it is 
very likely that precision experiments will discover 
lepton flavor violation via radiative loops. Dedic- 
ated LFV search experiments like the proposed 
/i — > eee experiment would then allow one to 
measure the LFV couplings of the new particles, 
complementary to the TeV scale experiments at 
the LHC. 

Conversely, in the case that no new physics (ex- 
cluding the SM Higgs boson 



114 



115 ) were dis- 



covered at the LHC, the discovery of LFV in preci- 
sion experiments is not excluded as e.g. rare muon 
decays are testing the mass scale > 1 PeV, three 
orders of magnitude higher than at LHC. 
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Decay Belle limit Babar limit Belle II proj. Belle II proj. SuperB proj. 1 



CHANNEL (Sab" 1 ) (50 atr 1 ) (75 ab" 1 ) 



r — > /i7 


4.5 


10" 8 


26 


4.4 


10" 8 


27 


10 


• io- 9 


|42 


43 


3 


io- 9 


42 


43 


1.8 • IO" 9 


96 




r — > e7 


12 • 


io- 8 




3.3 


io- 8 


27 


















2.3 • 10~ 9 


96 




T — > /i/i/i 


2.1 


io- 8 


34 


3.3 


io- 8 


28 


3 • 


IO" 9 


42 


43 


1 


io- 9 


42 


43 


2-10- 10 


96 




t — > eee 


2.7 


io- 8 


34 


2.9 


io- 8 


28 


















2-10- 10 


96 




t -4 /i?7 


2.3 


io- 8 


25 


15- 


io- 8 


33 


5- 


10~ 9 


42 


43 


2 


10 -9 


42 


43 


4-10- 10 


96 




r — > e?7 


4.4 


io- 8 


25 


16- 


io- 8 


33 


















6-10- 10 


96 






2.3 


io- 8 


35 


4.0 


io- 8 


31 


















2-10- 10 


96 




r -4 cfTg 


2.6 


IO" 8 


35 


3.3 


IO" 8 


31 


















2 • IO" 10 


96 





Table 3.3: Measured and projected limits on selected lepton flavour violating r decays (90% C.L.). 

1 The SuperB projections assumed a polarized electron beam; they also assumed that all backgrounds except initial state 
radiation can be suppressed to the desired level. The SuperB project was canceled in November 2012. 
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Chapter 4 



At 



The Decay h —> eee 



4.1 Kinematics 



The decay /i — > eee proceeds promptly. For dis- 
criminating signal and background, energy and 
momentum conservation can be exploited. The 
vectorial sum of all decay particle momenta should 
vanish: 

\ptot\ = = C 4 - 1 ) 

and the total energy has to be equal to the muon 
mass. 

The energies of the decay electrons (positrons) 
are in the range (0 — 53) MeV. All decay particles 
must lie in a plane and the decay is described 
by two independent variables in addition to three 
global rotation angles, which describe the orient- 
ation in space. 



4.2 Detector Acceptance 

The acceptance of the proposed /i — > eee exper- 
iment is determined by its geometrical accept- 
ance and energy coverage. For various coupling 
assumptions about the LFV amplitude, see also 
equation |2.1| the energy spectrum of the highest 
energy, E\, and lowest energy decay particles, 
, are shown in Figures 



E 



4.1 



and 4.2 



respect- 



ively. In order to achieve a high acceptance, the 
detector must be able to reconstruct tracks with 
momenta ranging from half the muon mass down 
to a few MeV with large solid angle coverage. 
The proposed experiment should cover the energy 
range > 10 MeV to provide acceptances of 50 % or 
more for all models. 




Figure 4.1: Energy distribution of the highest en- 



ergy positron in the decay fjr 



for dif- 



ferent effective LFV models. The solid red and the 
green lines correspond to pure four-fermion con- 
tact interaction models (no penguin) contribution. 



4.3 Backgrounds 



The final sensitivity of the proposed experiment 
depends on the ability to reduce backgrounds from 
various sources. Two categories of backgrounds 
are considered; irreducible backgrounds, such as 
fi + — > e + e + e^vv, which strongly depend on the 
granularity and resolution of the detector, and ac- 
cidental backgrounds that scale linearly or with 
the square of the beam intensity. 

In the following sections, the main background 
sources considered are discussed. 
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40 50 

E e min (MeV) 

Figure 4.2: Acceptance of the lowest energy de- 
cay electron (positron) for different effective LFV 
models as function of the minimum transverse mo- 
mentum. The solid red and green lines correspond 
to pure four-fermion contact interaction models 
(no penguin) contribution. 



4.3.1 Internal Conversions 

The decay fi — > eee vv o ccurs with a branching 
fraction of 3.4- 10" 5 ~ 



117 



. It can be distinguished 
from the /i — > eee process by making use of en- 
ergy and momentum conservation to reconstruct 
the undetected neutrinos; in order to separate the 
fi — > eee events from fj, — > eeevv events, the total 
momentum in the event is required to be zero and 
the energy equal to the muon rest energy. The 
branching fraction as a function of the energy cut 
of the /i —> eeevv process 
ure ■ 



4.3 Figures. 4.4 and 4.5 



116 is shown in Fig- 
show the energy spec- 
trum of all and the lowest energy electron from 
internal conversion decays, Figs. |4.6| and |4.7| the 
invariant masses of e + e~ combinations calculated 



with the matrix element from 116 . This process 



is the most serious background for the /i — > eee 
search and can only be resolved by a very good 
energy resolution. 

4.3.2 Michel Decays 

Using a beam of positive muons, one of the 
main processes contributing to accidental back- 
ground is that of the ordinary Michel decay fi + — > 
e + vv. This process does not produce a negat- 
ively charged particle (electron), which is one of 



10 
icr 1 

icr 1 
icr 1 

10" 1 
MO" 1 

icr 1 

10 2 
10- 



All internal conversion decays 



All decay electrons in detector acceptance ' 



99 100 101 102 103 104 105 

Visible Mass [MeV/c 2 ] 

Figure 4.3: Integrated branching fraction for the 
decay /i — > eeevv in dependence of the visible 
mass for all internal conversion decays and those 
with all three decay particles in the detector ac- 
ceptance. The matrix element was taken from 



116 . 




e*/e" Energy [MeV] 



Figure 4.4: Spectrum of electrons from internal 
conversion decays. 




Lowest energy of e*/e" [MeV] 



Figure 4.5: Spectrum of the electron with min- 
imum energy from internal conversion decays. 
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100 

JMeV] 



I 



Figure 4.6: Invariant masses of the two possible 
e + e~ combinations for internal conversion decays. 




Figure 4.7: Invariant masses of the two possible 
e + e~ combinations for internal conversion decays 
with a visible mass above 90MeV and the elec- 
trons and positrons in the detector acceptance 
(E > lOMeV, | cos (9 1 < 0.8). 



the main characteristics of the /i + — > e + e + e~ de- 
cay, and can therefore only contribute as potential 
background if a track is wrongly reconstructed. 
Other processes which "naturally" provide negat- 
ively charged tracks (electrons) are radiative de- 
cays with internal or external photon conversions 
or Bhabha scattering. 

4.3.3 Radiative Muon Decays 

The process fi + — > e +r ^vv (branching fraction 
1.4 • 10 -2 for photon energies above lOMeV 



117 



) 

can deliver an oppositely charged electron if the 
photon converts cither in the target region or in 



the detector. Contributions from conversions out- 
side of the target are greatly suppressed if a vertex 
constraint is applied and by minimizing the ma- 
terial in both the target and detector. Photon 
conversion in the target generates an event topo- 
logy similar to the radiative decay with internal 
conversion: [i — > eeeuv, which is discussed above. 

Due to the missing energy from the neutri- 
nos, this process mainly contributes to the acci- 
dental background in combination with an ordin- 
ary muon decay. 

4.3.4 Bhabha Scattering 

Positrons from the ordinary muon decay or beam- 
positrons can undergo Bhabha scattering with 
electrons in the target material, leading to an 
electron-positron pair from a common vertex. Due 
to the missing energy, this process mainly contrib- 
utes to the accidental background in combination 
with an ordinary muon decay. 

4.3.5 PlON DECAYS 

Certain decays of pions, especially 7r — > eeev 
(branching fraction 3.2 • 10 -9 |l!7|) a nd 7r — > fvyv 



(branching fraction 2.0 • 10 -4 117 ) with sub- 
sequent photon conversion arc indistinguishable 
from signal events if the momenta of the final 
state particles fit the muon mass hypothesis; a 
low pion contamination of the primary beam (es- 
timated to be in the order of 10~ 12 for the high 
intensity beamline), the small branching fraction 
and the small slice of the momentum is assumed 
to lead to negligible rates in the kinematic region 
of interest. 

4.3.6 Summary of Background Sources 

First simulation studies have been performed to 
calculate the different background contributions. 
Their results indicate that purely accidental back- 
grounds for ~ 10 9 muons stops per second are 
small for the proposed high resolution detector. 

The main concern are irreducible backgrounds, 
such as the process \i — > eeevv, which can only be 
reduced by a very good tracking resolution result- 
ing in total energy resolution of <je < 1 MeV for 
the aimed sensitivities < 10~ 15 . 
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The Mu3e Experiment 
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Chapter 5 




Requirements for Mu3e 



5.1 Goals of the Experiment 

The goal of the Mu3e experiment is to observe the 
process fi — > eee if its branching fraction is larger 
than 1CP 16 or otherwise to exclude a branching 
fraction of > 1CP 16 at the 90% certainty level. 
In order to achieve these goals, > 5.5 • 10 16 muon 
decays have to be observecj]] and any background 
mimicking the signal process has to be suppressed 
to below the 10 -16 level. The additional require- 
ment of achieving these goals within a reasonable 
measurement time of one year of data taking dic- 
tates a muon stopping rate of 2 • 10 9 Hz and a high 
geometrical acceptance and efficiency of the exper- 
iment. 

We plan to perform the experiment in two 
phases. The exploratory phase I will make use 
of existing muon beams at PSI and serve to com- 
mission the detectors, gain experience with the 
new technologies and validate the experimental 
concept, whilst at the same time producing a com- 
petitive measurement. The goal for this first phase 
is to reach a sensitivity of 10 -15 , thus pushing 
the existing limit by three orders of magnitude. 
For this level of sensitivity, the demands on the 
detector are somewhat relaxed, thus allowing for 
cross-checks between detectors also on analysis 
level or running without the full instrumentation. 
The lower data rates also will not require the full 
read-out and filter farm system. The second phase 
of the experiment on the other hand will aim for 
the ultimate sensitivity and thus require that the 
detector works as specified and a new beamlinc 
delivers > 2 • 10 9 Hz of muons. 

1 Assuming a total efficiency of 30%. 



The expected rate at an existing beamlinc is 
1 — 1.5 • 10 8 Hz of muons on target. In order to 
have a safety margin, we usually assume 2 • 10 8 Hz 
for phase I background studies, except where the 
running time is concerned. 

This proposal discusses the phase I experi- 
ment in detail and shows the path leading to full 
rate capability. We also discuss alternative ap- 
proaches. 

5.2 Challenges for the Experiment 

5.2.1 Backgrounds 

There are two kinds of backgrounds: Overlays of 
different processes producing three tracks resem- 
bling a (i — > eee decay (accidental background) 
and radiative decays with internal conversion (in- 
ternal conversion background) with a small energy 
fraction carried away by the neutrinos. Accidental 
background has to be suppressed via vertexing, 
timing and momentum measurement, whereas 
momentum measurement is the only handle on in- 
ternal conversion. 

5.2.2 Geometric acceptance 

For a three-body decay with a priori unknown kin- 
ematics such as /i — y eee, the acceptance has to 
be as high as possible in order to test new phys- 
ics in all regions of phase space. There are two 
kinds of acceptance losses, losses of tracks down- 
stream or upstream, where beam entry and exit 
prevent instrumentation, and losses of low trans- 
verse momentum tracks, which do not transverse 



22 




An Experiment to Search for the Decay fi — > eee 



.2 1(r « 






1 10" 2 












§,10 13 






I 10 14 






i io- 15 




— -^"^ 5 sigma around signal 


I itr 16 




4 sigma around signal 


| 10"" 




3 sigma around signal 


o 10 




2 sigma around signal 


t 10 19 




1 sigma around signal 


i 10 -20 


IT i i i i i i i i i i 





0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

Reconstructed Mass Resolution [MeV/c 2 ] 



Figure 5.1: Contamination of the signal region 
(one sided cut) with internal conversion events as 
a function of momentum sum resolution. 

a sufficient number of detector planes, and are not 
reconstructed. 

5.2.3 Rate capability 

The Mu3e detector should be capable of running 
with 2 • 10 9 Hz of muon decays. This poses chal- 
lenges for the detectors, the data acquisition and 
the readout. 



5.2.4 Momentum resolution 

The momentum resolution directly determines to 
what level internal conversion background can be 
suppressed and thus to which level the experiment 
can be ran background free. In order to reach a 
sensitivity of 1CP 16 with a 2er cut on the recon- 
structed muon mass, the average momentum res- 
olution has to be better than 0.5 MeV. For the 
phase I experiment aiming at 10 -15 , this require- 
ment is relaxed to 0.7 MeV, see Figure [5~T| 

5.2.5 Vertex resolution 

Keeping apart vertices from different muon decays 
is a key tool in suppressing accidental background. 
The vertex resolution is essentially determined by 
the amount of multiple scattering (and thus ma- 
terial) in the innermost detector layer. Ideally 
the vertex resolution is sufficient to eliminate al- 
most all combinatorial backgrounds; for the phase 
I rates, this appears achievable, whereas in the 
phase II experiment, very good timing is needed 
in addition. 

5.2.6 Timing resolution 

Good timing is essential for reducing combinator- 
ial background at rates which lead to more than 
about 10 muon decays per frame on average. 
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Chapter 6 




Experimental Concept 



The Mu3e detector is aimed at the background 
free measurement or exclusion of the decay /i — > 
eee at the level of 10 -16 . As discussed in more 
detail in the preceding chapter [5] these goals re- 
quire to run at high muon decay rates, an excel- 
lent momentum resolution in order to suppress 
background from the internal conversion decay 
fi — > eeevv and good vertex and timing resolu- 
tion in order to efficiently suppress combinatorial 
background. 

We intend to measure the momenta of the muon 
decay electrons in a solenoidal magnetic field us- 
ing a silicon pixel tracker. At the electron ener- 
gies of interest, multiple Coulomb scattering in de- 
tector material is the dominating factor affecting 
momentum resolution. Minimizing this material 
in the active detector parts is thus of utmost im- 
portance. 



The proposed detector consists of an ultra thin 
silicon pixel tracker, made possible by the High- 
Voltage Monolithic Active Pixel (HV-MAPS) 
technology (see chapter 10 1. Just four radial lay- 



ers around a fixed target in a solenoidal magnetic 
field allow for precise momentum and vertex de- 
termination. Two timing detector systems guar- 
antee good combinatorial background suppression 
and high rate capabilities. 

The Mu3e experiment is designed to have a 
sensitivity four orders of magnitude better than 
the current limit on /i — !> eee (1CP 12 ), so it is reas- 
onable to plan for a staged detector design, with 
each stage roughly corresponding to an order of 
magnitude improvement. 
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6.1 Momentum Measurement with 
Recurlers 

Due to the low momenta of the electrons from 
muon decay, multiple scattering is the dominat- 
ing effect on momentum measurement. With our 
fine-grained pixel detector, we are thus in a regime 
where scattering effects dominate over sensor res- 



olution effects, see Figs. |6.f| and |6.2| Thus adding 
additional measurement points does not necessar- 
ily improve the precision. 

The precision of a momentum measurement de- 
pends on the amount of track deflection fl in the 
magnetic field B and the multiple scattering angle 
©MS, see Figure [6~3| to first order: 



9 



It is 



P 



o 



(6.1) 



So in order to have a high momentum precision, 
a large lever arm is needed. This can be achieved 
by moving tracking stations to large radii, which 
however compromises the acceptance for low mo- 
mentum particles. In the case of muon decays, all 
track momenta are below 53MeV and all tracks 
will thus curl back towards the magnet axis if 
the magnet bore is sufficiently large. After ex- 
actly half a turn, effects of multiple scattering on 
momentum measurement cancel in first order, see 
Figure |6.4| To exploit this feature we optimized 
the experimental design specifically for the meas- 
urement of re-curling tracks, leading to a narrow, 
long tube layout. 

Measuring the momentum from bending outside 
of the tracker also allows us to place timing de- 
tectors inside, without strongly affecting the res- 
olution. 

6.2 Baseline Design 

The proposed Mu3e detector is based on two 
double layers of HV-MAPS around a hollow 
double cone target, see Figures [6~5] and [676] The 
outer two pixel sensor layers are extended up- 
stream and downstream to provide precise mo- 
mentum measurements in an extended region with 
the help of re-curling electrons. The silicon de- 
tector layers (described in detail in chapter [To| 
arc supplemented by two timing systems, a scin- 
tillating fibre tracker in the central part (see 
chapter 111 and scintillating tiles (chapter 12 1 in- 
side the recurl layers. Precise timing of all tracks 




Figure 6.3: Multiple scattering as seen in the plane 
transverse to the magnetic field direction. 




Figure 6.4: Multiple scattering for a semi-circular 
trajectory. 
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is necessary for event building and to suppress ac- 
cidental combinatorial background. 

The entire detector is built in a cylindrical 
shape around a beam pipe, with a total length of 
approximately 2 m, inside a 1 T solenoid magnet 
with lm inside diameter and 2.5 m total length 
(chapter [8]) . In the longitudinal direction the de- 
tector is sub-divided into five stations, the central 
detector with target, inner silicon double layer, 
fibre tracker and outer silicon double layer, and 
two forward and backward recurl stations with 
two silicon recurl layers surrounding a tile tim- 
ing detector. In order to separate tracks coming 
from different muon decays, the target has a large 
surface with 10 cm length and 2 cm diameter. The 
target shape is a hollow double cone, see chapter 
[9] Around the target the two inner silicon pixel 
layers, also referred to as the vertex layers, cover 
a length of 12cm. The innermost layer will have 
12, the second one 18 sides of 1cm each, corres- 
ponding to an average radius of 1.9 cm and 2.9 cm, 
respectively. The inner silicon layers are suppor- 
ted by two half cylinders made from 25 pm thin 
Kapton foil mounted on plastic end pieces. All 
silicon sensors arc thinned to 50 um, resulting in a 
radiation length of X/Xo < 0.1% per layer. The 
detector will be cooled with gaseous helium. 

The hit information from the silicon sensors is 
read out at a rate of 20 MHz using timestamps 
providing a time resolution of 20 ns. 

The fibre tracker sits inside silicon pixel layer 
three at around 6 cm, providing timing informa- 
tion for decay positrons and electrons. It is com- 
posed from three to five layers of 250 pm thick 
36 cm long scintillating fibres, see Figure [6~6l The 
fibre tracker is read out by fast silicon photo mul- 
tipliers and can provide timing information with 
< 1 ns accuracy. 

The silicon pixel layers three and four are just 
outside the fibre tracker at a mean radius of 7.6 cm 
and 8.9 cm. The active area has a cylindrical 
shape of 36 cm length. The layer three has 24 
sides, layer four 28 sides of 1.9 cm width each. 
Both outer layers are constructed as modules of 4 
sides, six modules for layer three and seven mod- 
ules for layer four. Similar to the inner two layers 
the mechanical frames of these modules are build 
from 25 pm Kapton foil with plastic end pieces. 

Copies of silicon pixel layer three and four are 
also used in the recurl stations. Two recurl sta- 
tions each are covering the upstream and down- 
stream regions. These recurl stations add fur- 



ther precision to the momentum measurement of 
the electrons, see section |6.1| While the silicon 
layer design is (almost) identical to the central 
part, the timing detector in the recurl region can 
be much thicker compared to the fibre tracker, 
as the particles can and should be stopped here. 
This is done by using scintillating tiles of about 
7.5 x 7.5 x 5 mm 3 size. These tiles provide a much 
better time resolution than the thin fibre tracker 
in the center. Following the dimensions of the 
recurl silicon layers, the tile station have a act- 
ive length of 36 cm and a cylindrical shape with 
a radius of ss 6 cm. All central detector com- 
ponents are mounted on spokes providing a light 
stiff support. The recurl silicon layers and tiles 
are mounted on the beam pipe support. In addi- 
tion to the silicon and scintillating tile sensors the 
beam pipe support also carries the services and 
the PCBs equipped with the front-end electronics 



(chapter 131. Signal and power connection to the 



silicon layers is provided by flex prints which are 
also part of the mechanical support of the silicon 
sensors. 



6.3 Building up the Experiment 

One of the advantages of the design concept 
presented is its modularity. Even with a par- 
tial detector, physics runs can be taken. The 
full instrumentation is only required for achieving 
the final sensitivity of 10~ 16 at muon rates above 
1 • 10 9 Hz. On the other hand, in an early commis- 
sioning phase at smaller muon stopping rates, the 
detector could run with the central silicon detector 



only (see Figure 6.7 1. The silicon detectors of the 



recurl stations are essentially copies of the central 
outer silicon detector; after a successful commis- 
sioning of the latter, they can be produced and 
added to the experiment as they become avail- 
able together with the connected tile detectors. 
The fibre tracker can also be added later, since 
it is only needed to resolve combinatorial back- 
ground at higher event rates and track multiplicit- 
ies. The loss of momentum resolution due to mul- 
tiple scattering at the additional material of the 
fibre tracker will be fully compensated by the im- 
proved momentum measurement with re-curlers. 
The configuration with two recurl stations (Fig- 
defines a medium-size setup, well suited 



6. 



for phase I running. The configuration with four 
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Figure 6.5: Schematic view of the experiment cut along the beam axis in the phase II configuration. 




Figure 6.6: Schematic view of the experiment cut transverse to the beam axis. Note that the fibres are 
not drawn to scale. 



recurl stations (Figure 6.9 1 defines the full setup 
for phase II running. 



In the following sections, the experimental con- 
figurations for running at the existing 7rE5 beam- 
line (the Phase I Experiment) and the final de- 
tector for running at > 1 • 10 9 Hz muon stopping 
rate (the Phase II Experiment) are outlined. 



6.4 The Phase I Experiment 

The phase I of the Mu3e experiment will start with 
a minimum configuration (phase IA detector) with 
the target regions surrounded by double layers of 
inner and outer silicon pixel detectors, see Figure 
|6.7| This configuration defines the minimal config- 
uration as it allows to determine the momentum, 
the vertex position and the time of the decay pre- 
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Figure 6.7: Minimum detector cofiguration for early commissioning with central silicon only (phase IA). 




Figure 6.8: Detector with one set of recurl stations for physics runs and tile detector commissioning 
(phase IB). 




Figure 6.9: Final detector with two sets of recurl stations for high intensity physics runs (phase II). 
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cise enough to produce very competitive physics 
results with a sensitivity down to O(10" 14 ). It is 
foreseen to run in the first year in this configura- 
tion at a muon stopping rate on target at around 
2 • 10 7 Hz. The number of decays in one readout 
frame of the pixel tracker of 50 ns will be around 
one on average and combinatorial background can 
be suppressed with the help of the vertex recon- 
struction. The precision of the momentum res- 
olution will be somewhat limited, as most tracks 
do not recurl within the instrumented volume, see 



chapter 17 



In the phase IB the detector will be complemen- 
ted by the first pair of recurl stations, the corres- 
ponding tile detectors and the fibre tracker, see 
Figure |6.8| Adding the recurl stations will sig- 
nificantly enhance the momentum resolution and 
thus improve the suppression of internal conver- 
sion background. The insertion of the fibre tracker 
and the tile detector stations gives a much better 
time resolution in comparison to the silicon pixel 
only. The fibre tracker will deliver a time resol- 
ution of about 200-300 ps, while the tile detector 
will have < 100 ps resolution for the tracks passing 
the recurl stations. The high time resolution will 



allow running at the highest possible rate at the 
7rE5 muon beam line at PSI of « 1 • 10 8 Hz. The 
sensitivity reach in this phase of the experiment of 
0(1O~ 15 ) will be limited by statistics only, limited 
by the available muon decay rate. 

6.5 The Phase II Experiment 



A new high intensity muon beam line |118| deliver- 
ing w 2 • 10 9 Hz muon stops is crucial for the phase 
II of the proposed experiment. To fully exploit the 
new beam facility the limited detector acceptance 
at phase IB will be further enhanced by adding 
another a second pair of recurl and tile detector 
stations, see Figure |6.9| These extra stations will 
allow to measure precisely the momentum of all 
particles in the acceptance of the inner tracking 
detector. At the same time the extra tile detector 
stations with their high time resolution and small 
occupancy will help to fight the increased com- 
binatorics at very high decay rates. The com- 
bined performance of the final detector setup to- 
gether with the high stopping rate will allow to 
search for the fj, — > eee decay with a sensitivity of 
B(/i -> eee)< 10' 
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Chapter 7 




Muon Beam 



7.1 General Beam Requirements 

The general beam requirements for a high intens- 
ity, low-energy, stopped muon coincidence exper- 
iment such as Mu3e are six-fold: an abundant 
supply of low-energy surface muons (from stopped 
pion decay at rest, at the surface of the produc- 
tion target 119] ) capable of being stopped in a 
thin target; high transmission optics at 28MeV/c, 
close to the kinematic-edge of stopped pion de- 
cay and hence close to the maximum production 
rate of such muons, as shown in Figure |7.1| a 
small beam emittance to minimize the stopping- 
target diameter; a momentum-byte of less than 



10% with an achromatic final focus, allowing an 
almost monochromatic beam with a high stopping 
density, to be stopped in a minimally thick tar- 
get; minimization and separation of beam-related 
backgrounds such as beam e + originating from ir°- 
decay in the production target, or decay particles 
produced along the beam line and finally minim- 
ization of material interactions in the beam, for 
example such as those in windows, thus requir- 
ing vacuum or helium environments to keep the 
multiple scattering under control. 



7.2 Beam for phase I running 



4000 



Z- Branch Momentum Spectrum 




30 31 32 33 
Muon Momentum [MeV/c] 

Figure 7.1: 7rE5 measured muon momentum spec- 
trum with fully open slits. Each point is obtained 
by optimizing the whole beam line for the cor- 
responding central momentum and measuring the 
full beam-spot intensity. The red-line is a fit to 
the data, based on a theoretical p 3 5 behaviour, 
folded with a Gaussian resolution function corres- 
ponding to the momentum-byte plus a constant 
cloud- muon background. 



As previously outlined, a multi-staged detector 
configuration will be sought for phase I running, 
this in turn requires muon beam intensities ran- 
ging between (10 7 -10 8 ) muons/s for the initial 
phase I A (central detectors), while for phase IB 
a maximal intensity close to 1 • 10 8 muons/s will 
be sought. The quoted maximum beam intens- 
ity which includes a phase space reduction factor 
due to collimation in the central region of the 
Mu3 magnet is based on measured intensities at 
the centre of the MEG detector, without a de- 
grader and normalized to a proton beam current 
of 2.3 mA. This demand for the highest intensities 
necessitates the selection of only one possible facil- 
ity in the world, namely the 7rE5 channel at PSI. 
Based on the experience gained in the design of the 
MEG beam line, a similarly developed concept is 
also envisaged for phase I of Mu3e. This should 
allow the required muon intensity to be achieved. 
However, since the area is likely to be shared with 
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Mu3e Phase I Schematic 




"Z"-Channel 



Figure 7.2: Mu3e potential beam line layout in the front-part of the 7rE5 area. 



Separator Scan 




-500 500 1000 1500 2000 2500 



Separator DAC value 

Figure 7.3: Separator scan plot, measured post 
collimator system. The black dots represent on- 
axis, low-threshold intensity measurements at a 
given separator magnetic field value during the 
scan, for a fixed electric field value of — 195 kV. 
The red curve represents a double Gaussian fit to 
the data points, with a constant background. A 
separation of 8.1 muon beam is found, corres- 
ponding to more than 12 cm separation of beam- 
spot centres at the collimator. 



other experiments, that of MEG (R-99-05), and 
the Lamb-shift experiment (R-98-03) a compact 
muon beam line designed specifically to fit into 
the front-part of the 7rE5 area is under design. 
This would not only allow the beam line elements, 
such as the Wien-filter, triplet I and II, plus the 
MEG collimator system to be used by this exper- 
iment, but would also allow access to the MEG 
detector during running periods, by means of pla- 
cing a shielding wall just upstream of the MEG 
detector, as previously adopted during the run 
periods of experiment R98-03. 

Figure [73] shows the potential area layout adop- 
ted. Surface muons of 28MeV/c will be extracted 
from the 7rE5 "Z-channel" and the initial part of 
the current MEG beam line, including: Triplet 
I, the 200 kV crossed-field Wien-filter, Triplet II 
and the collimator system. This combination of 
elements allows for an optimal beam correlated 
background suppression, as demonstrated in Fig- 
ure |7.3| which shows the separation quality post 
collimator, between muons and beam positrons for 
the above mentioned section of beam line. Due to 
the severe restrictions imposed by space, a match- 
ing section, including two dipole magnets of 90° 
and 66° bending angles respectively, with an in- 
termediate quadrupole doublet or triplet, is en- 
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visaged. A final doublet /triplet or intermediate 
transport solenoid would in turn couple the beam 
line to the Mu3e superconducting magnet. The 
beam line vacuum is currently planned to end 
close to the centre of the target and includes col- 
limation to match the beam-spot to the target size 
and prevent beam interactions from occurring dir- 
ectly in the small radii inner silicon layers. The 
expected usable muon intensity at the Mu3e tar- 
get is between 0.7-1.0 • 10 8 muons/s. 

Although it is understood that simultaneous 
running of prospective 7rE5 experiments is not 
possible, it is nevertheless clear that the Mu3e 
phase I beam line will in fact benefit from the 
availability of MEG beam elements upstream of 
the detector and that this option, together with 
provision of the available standard PSI magnets 
currently in storage or potentially sharable, would 
cover most of the beam line requirements, except 
that of a short coupling solenoid, for which a po- 
tential solution is also currently under study and 
possible dipole vacuum chamber modifications. 

7.3 High intensity muon beamline 
for phase II running 

In order to reach the ultimate sensitivity goal 
of 0(1O -16 ) for the phase II experiment, an im- 
pulsed muon stopping rate in the GHz region is 
required. As demonstrated in Table 7.1 there 



are no such (pulsed or unpulsed) high-intensity 
sources of muons currently available world-wide. 
Future intensity frontier facilities are however in 
the planning in the US and Japan and are also 
associated with LFV-experiments, more specific- 
ally Mu2e and Project X in the U.S. (461P] and 



the COMET and PRIME/PRISM experiments in 
Japan 45 101] , To meet the needs of such ex- 



periments a totally new concept is therefore ne- 
cessary. One such concept, which is still in its 
infancy, though is proving to be a promising can- 
didate, is the HiMB project at PSI [118] , a next- 
generation high-intensity muon beam, currently 
under study. A detailed feasibility study is due to 
start at PSI in 2013. This concept would provide 
the basis for a new Mu3e beam line for the phase II 
measurements, based on the production of surface 
muons from the Swiss Spallation Neutron Source's 
(SINQ) spallation target window. 



with the proton beam injected from below is 
shown in Figure |7.4| The characteristics of the 
source, which resembles a medium-flux reactor, 
are that the protons are injected from below and 
defocussed onto a double-layered aluminium win- 
dow separated by a D 2 cooling layer, before be- 
ing stopped in the target, a "cannelloni" construct 
of lead-filled zircaloy tubes. 

The HiMB project plans to extract the down- 
ward travelling muons produced in the aluminium 
window via a two-stage channel, the first stage, 
a solenoidal one, which uses the same beam- 
tube as the upward travelling protons and ex- 
tracts the muons, in the opposite direction, to a 
large collection solenoid connecting to the second 
stage, a conventional dipole and quadrupolc chan- 
nel planned for the empty service cellar under the 
SINQ target. The general layout of both the pro- 
ton channel and the service area are shown in Fig- 
ure [723 

There are several advantages of this concept, 
which would lead to a substantial enhancement 
compared to target E, namely: the increased num- 
ber of primary proton interactions since 70 % of 
the beam stops in the target; a much larger pion 
energy range of up to 150 MeV can be exploited 
in the case of SINQ, above which the high en- 
ergy tail of the pion production cross-sections be- 
comes negligible, in the case of Target E this limit 



is around 45 MeV 121 122 ; a substantially larger 



The layout of the source in the SINQ hall 120 



together with a schematic diagram of the source 



pion-production volume contribution compared to 
Target E and finally a significantly larger surface 
muon production volume. 

Realistic Monte-Carlo studies were undertaken 
together with M. Wohlmuther (Head of the Target 
Development Group at PSI) using the Los Alamos 
Laboratory MCNPX code (Monte Carlo n-particle 
extended code), used also to design the SINQ tar- 
get. A simulation of the surface muon produc- 
tion rate was made using the complete model of 
the SINQ target environment. Based on a total 
of 4 • 10 8 generated upward moving protons, cor- 
responding to the measured 2D beam profile at 
SINQ, a complete particle tracking was done us- 
ing three different event generators. For surface 
muons, the simulated fluences were determined for 
the conditions of a particle leaving the target with 
the correct energy, travelling downwards within 
the beam-pipe and crossing a horizontal plane 
25 cm below the window. This is shown schemat- 
ically in Figure |7.6| The calculated fluences from 
the three event generators agreed to within 35% 
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Figure 7.4: (Right) The SINQ spallation neutron source as seen in the experimental hall. The tower 
in the foreground, shown in section with the proton beam incident from below, (left) contains the 
main neutronics components, the spallation source, a Pb/Zr/Al structure cooled using D2O, the D2O 
moderator tank and reflector shield and a cold-source of solid deuterium. The upper-part of the tower 
deals with the cryogenic connections of the source. 




Figure 7.5: Shows the layout of the proton beam line at SINQ as well as showing the possibility for 
extracting the muons into the empty cellar region below the SINQ target. The blue elements are dipolc 
bending magnets and the red elements focussing quadrupole magnets. Since the proton and muons are 
like-sign charged particles but travelling in opposite directions they will be bent in opposite directions, 
allowing the muons to access the empty cellar region shown. This is where a muon beam line could be 
placed. 
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Laboratory/ Energy/ Present Surface Future estimated 

Beam line Power fi + rate (Hz) /i + / fjb~ rate (Hz) 

PSI (CH) (590 McV, 1.3 MW, DC) 

LEMS " 4 • 10 s 

7rE5 " 1.6 ■ 10 8 

HiMB (590 MeV, 1 MW, DC) 4 ■ 10 10 (/i+) 



J-PARC (JP) (3 GeV, 1 MW, Pulsed) 

currently 210 KW 

MUSE D-linc " 3 • 10 7 

MUSE U-line " 2 ■ 10 s (At+) (2012) 

COMET (8 GeV, 56 kW, Pulsed) lO 11 ^") (2019/20) 

PRIME/PRISM (8 GeV, 300 kW, Pulsed) lO 11 " 12 ^") (> 2020) 



FNAL (USA) 

Mu2e 

Project X Mu2e 


(8 GeV, 25 kW, Pulsed) 
(3 GeV, 750 kW, Pulsed) 




5 • 10 10 (/O (2019/20) 
2 ■ 10 12 (/i-) (> 2022) 


TRIUMF (CA) 

M20 


(500 MeV, 75 kW, DC) 


2 ■ 10 6 




KEK (JP) 

Dai Omega 


(500 MeV, 2.5 kW, Pulsed) 


4- 10 5 




RAL -ISIS (UK) 

RIKEN-RAL 


(800 MeV, 160 kW, Pulsed) 


1.5 ■ 10 6 




RCNP Osaka Univ. (JP) 

MUSIC 


(400 MeV, 400 W, Pulsed) 
currently max 4W 




10 8 (At+) (2012) 
means > 10 11 per MW 


DUBNA (RU) 

Phasatron Ch:I-III 


(660 MeV, 1.65 kW, Pulsed) 


3 ■ 10 4 





Table 7.1: Currently running muon beam facilities around the world used for particle physics experiments 
and materials science /iSR investigations. Also shown are the planned next-generation facilities designed 
for cLFV experiments, together with an estimate of the starting date. The PSI HiMB solution is currently 
only under study and is included purely for completeness. 



of each other. Based on the standard event gener- 
ator, which also has the smallest statistical uncer- 
tainty, a summed fluence (E^ 4.12 MeV) of surface 
and sub-surface muons of 1 ■ 10 11 muons/s is cal- 
culated at a proton current of 3 mA on target E, 
which corresponds to 2.1mA on SINQ. However, 
on the assumption that the proton current on Tar- 
get E will only rise to a maximum of 2.4mA in 
the future, a value that has already been achieved 
during routine test periods since 2010 and taking 
into account the variation of event generators in 
the simulation, a conservative estimated fluence 
of (3 ± 1) ■ 10 10 good surface muons within a 10 % 



FWHM momentum-byte could clearly be extrac- 



ted from the SINQ target (c.f. Table |7.2|. 

Conclusions seen from the source point of view 
look very promising and providing the beam can 
indeed be transmitted without significant losses, 
the phase II rates of 2 GHz muons could be real- 
ized. Finally, the current extraction principle, 
to be studied extensively in the feasibility study, 
is demonstrated in Figure 7.7 Muons originat- 



ing from the target window are to be guided in 
a downward direction using a low-field guiding 
solenoid. Since the muon momentum is about a 
factor of 42 times smaller than that of the protons, 
the low-field for the muons should not dramat- 
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Figure 7.6: The Monte-Carlo model for the SINQ 
target with all components taken into account. 
Light grey is the aluminium window, with the 
D2O cooling channel in the middle. Protons are 
generated according to a measured 2-D distribu- 
tion at the dark red plane. Surface muons leav- 
ing the target at the correct energy and cross- 
ing the light-green plane, 25 cm below the window 
and moving downwards within the beam-pipe, are 
counted. 



Pa 


AP/P 


Rate 


MeV/c 


% FWHM 


Hz 


28 


Full 


(7±1)-10 10 


28 


10 


(3±1)-10 10 


26 


10 


(3±1)-10 10 



Table 7.2: Estimated surface and sub-surface 
muon rates based on a proton current of 2.4mA 
on Target E and full transmission efficiency 
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Figure 7.7: Schematic of HiMB muon extraction 
principle, with a guiding solenoid, followed by a fo- 
cussing solenoid to satisfy proton and muon trans- 
mission. The extraction is done in the fringe-field 
of the AHO magnet, with a strong- focussing col- 
lection solenoid. This is followed by a conventional 
beam line consisting of a dipole magnet and quad- 
rupole channel. 



ically affect the protons. The present radiation- 
hardened, defocusing quadrupolcs QTH 31 and 
32, just below the SINQ target, must be replaced 
by a rotationally symmetric clement such as a 
solenoid, since the muons will not traverse a set 
of quadrupoles which only focus alternately in 
the horizontal and vertical planes but have field 
strengths that are 42 times too high for the muons. 
The extraction will be done in the fringe-field of 
the last dipole magnet AHO, which means that 
the following strong-focussing collection solenoid 
must be placed close to the "fan-coupling" in or- 
der to fully collect the beam. Following this solen- 
oid is a dipole magnet whose function is to bend 
the beam onto the horizontal plane, where a con- 
ventional secondary beam line quadrupole chan- 
nel could be constructed. As the current cellar 
ends within just a few meters of the above SINQ 
hall wall, one could envisage bending the beam 
upwards again at the end of the cellar and ex- 
tracting to a hall, exterior to SINQ, on the East- 
side. This would imply a relatively long beam line 
of order (30-35) m, which is not too problematic 
from the muon loss (decay-in-flight) point-of-view 
where 80% transmission for 28MeV/c muons is 
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Figure 7.8: Survival probability for muons at the kinematic edge of 29.79 MeV/c (green curve ), surface 
muons of 28MeV/c (blue curve) and sub-surface muons of 25MeV/c (red curve). Also shown (orange 
curve) is the survival probability for pion contamination at 28MeV/c, which is at the 1CP 12 level. 



expected at 40 m, as shown in Figure |7.8| Also 
shown arc similar plots for the maximum mo- 
mentum at the kinematic edge of pion-decay, as 
well as for a lower momentum sub-surface muon 
beam. Finally, the beam would naturally be free 
of pion contamination to a level of about 10 -12 , 
which is also important for backgrounds relevant 
to the (i eee decay such as ir — > eeev and 

The HiMB project is in its infancy at present 
and there are many aspects that will be studied 
in detail within the scope of the planned feasib- 
ility study to show the feasibility of this next 



generation high-intensity beam line. The initial 
step taken concerning the muon source intensity, 
the basis for the HiMB feasibility study, has been 
shown to be very promising. The next major steps 
to be studied are the optical extraction of the 
muons from the proton beam and the solenoidal 
replacement of the final proton defocussing quad- 
rupole doublet in front of the SINQ target, while 
still maintaining the safety restraints on the SINQ 
target. The implementation of such a concept into 
the SINQ environment could only coincide with a 
major SINQ shutdown, which is currently planned 
for the period of 2016-2017. 
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Magnet 



The magnet for the Mu3e experiment has to 
provide a homogeneous solenoidal magnetic field 
for the precise momentum determination of the 
muon decay products. In addition it will also serve 
as beam optical element guiding the muon beam 
to the target. The basic parameters of the super- 
conducting solenoid magnet, which are currently 
being specified for the preparation of order place- 
ments, are given in Table |8.1| The outer dimen- 
sions include also an iron field shield. 

The nominal magnetic field strength is 1 T in 
the central part, providing the optimum bending 
radius in terms of resolution for the proposed ex- 
perimental design. A higher magnetic field would 
lead to a loss of acceptance as the low momentum 
particles would not reach the central outer pixel 
layers (see Figure 8.1 1. A lower magnetic field 



would lead to less magnetic deflection at constant 
multiple scattering, leading to worse momentum 
resolution (see Figure 8.2 1. For systematic studies 



and to allow for possible reuses of the magnet for 
other experimental measurements, the field can be 
varied between 0.8 and 2T. 

The dimensions of the cylindrical warm bore 
of the magnet are 1 m in diameter and > 2 m 
in length. The minimum diameter is given by 
four times the bending radius of the highest mo- 
mentum (53MeV/c) decay products at the lowest 
possible field of B = 0.8 T plus the target dia- 
meter. In addition the detector support and ex- 
traction rail system has to be taken into account 
when choosing the warm bore diameter. 

The total length is a compromise between geo- 
metric acceptance for recurling particles and the 
very tight space constraints for the phase I exper- 



imental area at 7rE5. In principle a longer solen- 
oidal magnet would provide an intrinsically more 
homogenous field. At both ends of the magnet it is 
foreseen to have full access by means of removable 
flanges. 

While the ideal magnet would have a constant 
field throughout the inner volume, real solenoid 
magnets show a drop in field to 50 % at the end 
of the coil. The simplest solution would be a 
longer magnet, which however does not fit inside 
the phase I area. Another possibility is to intro- 
duce correction coils at both ends of the magnet, 
such that the high field region can be extended. 
The insertion of several compensating coils would 
make the magnet system more complex both in 
construction and operation due to the need of ad- 
ditional current settings and power supplies. At 
present the baseline magnet concept foresees three 
equal coils with a single power supply. The field 
change along the z-axis has to be taken into ac- 
count for the reconstruction of tracks in the recurl 



Magnet parameter 



Value 



field for experiment 

field range 

warm bore diameter 

warm bore length 

field description AB/B 

field stability AB/B (100 days) 

outer dimensions: length 
width 
height 



IT 

0.8 < 1 < 2T 
1 m 

> 2m 

< 10- 4 

< 10~ 4 

< 2.5 m 

< 2.5 m 

< 3.5 m 



Table 8.1: Properties of the Mu3e magnet 
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Figure 8.1: Reconstruction efficiency as a function 
of track transverse momentum for different mag- 
netic fields, with recurl stations and without fibre 
detector. 
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Figure 8.2: Momentum resolution as a function 
of track momentum for different magnetic fields, 
with recurl stations and without fibre detector. 



stations by using a look up table for the field map 
plus interpolation between these points. Choos- 
ing the right granularity for the look-up table a 
linear interpolation of the field will be enough to 
reach an approximation of AB/B < 10~ 4 . For 
the fast online selection of events the assumption 
of a constant field in the active part of the exper- 
iment will be sufficient. Though the assumption 



of a constant (maximum) field leads to an system- 
atic bias towards larger momenta and an increase 
of online selected background events from internal 
radiative muon decays with internal conversions, 
no signal events would be lost. 

The superconducting magnet is made from 
three coils of equal size, which has advantages over 
one long coil in terms of mechanical stability. The 
small dips in the magnetic field can be treated 
numerically in the same way as the roll-off of the 
field to the ends of the magnet. The choice for 
the superconducting wires or conductors will be 
driven by commercial availability, since standard 
components allow for the desired 2 T maximum 
field strength. A warm normal conducting mag- 
net is no option because of size, cost (copper price) 
and operational stability. Superconducting mag- 
nets have an intrinsic immunity against absolute 
field changes, as they have to run at a constant 
(low) temperature. If feasible in terms of number 
of cooling compressors, a dry cooled system will 
be chosen. 

There will be a magnetic shielding around the 
magnet. The shielding is required since the experi- 
mental hall is densely populated with other exper- 
iments and infrastructure. Also for the read-out of 
the proposed experiment it will be much easier to 
work in a low field environment. A beneficial side 
effect of the shielding is a gain of field homogen- 
eity inside the magnet and less field dependence 
on variation of outside parameters. 

The long term stability of the magnetic field 
should be AB/B < 10" 4 over each 100 day data 
taking period. This can be achieved by using state 
of the art magnet power supplies and by perman- 
ent measurement of the absolute field with a hall 
probe inside the experiment. 

The cool-down time for a system of the projec- 
ted size will be one week and the ramp time will be 
in the order of one hour. The number and power of 
the dry compressors will be chosen to fulfill these 
requirements, in the case of a dry cooled magnet. 

The DO Magnet |123| fulfills most requirements 
of the future Mu3e magnet and serves as a proto- 
type for the magnet design process. 
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Chapter 9 




Stopping Target 



The main challenge for the design of the stop- 
ping target is to optimize the stopping power on 
one hand and to minimize the impact on the track 
measurement on the other hand. Therefore the 
stopping target should contain just enough ma- 
terial in the beam direction to stop most of the 
29MeV surface muons but should be as thin as 
possible in the flight direction of decay electrons 
measrued in the detector acceptance. Usage of a 
low Z material is advantageous as tails from large 
angle Coulomb scatterering are suppressed. In ad- 
dition, the decay vertices should be spread out as 
widely as possible in order to reduce accidental co- 
incidences of track vertices and to produce a more 
or less even occupancy in the innermost detector 
layer. 



9.1 Baseline Aluminium Design 




30 Mm A! 



Al 




A = 3204 mm 2 



Figure 9.1: Dimensions of the baseline design tar- 
get. Note that the material thickness is not to 
scale. 



These requirements can be met by a hollow 
double cone target a la SINDRUM |15l|94l. In our 



baseline design (see Figure 9.1 ), the target is made 
from 30 um of aluminium in the front part and 
80pm aluminium in the back part, with a total 
length of 100 mm and a radius of 10 mm. This res- 
ults in an total area of 3204 mm 2 and an effective 
target thickness in beam direction of 560 um cor- 
responding to 0.063 radiation lengths X of Alu- 
minium. The target can be suspended from the 
innermost tracking layer by e.g. nylon fishing wire 
(which we assume in the simulation) and which 
does not significantly add material in the beam 
line. 



In the Geant4 |124[ simulation (see 15), about 
83.3% of the muon^j that reach the target are 
stopped. Obviously, this fraction can be increased 
by adding material, which will however lead to 
additional multiple scattering and thus a reduced 
momentum resolution. For the phase I experi- 
ment, where muon rates rather than momentum 
resolution is limiting the sensitivity, a thicker tar- 
get could be envisaged. 

Stopping 2 • 10 9 Hz muons in the target corres- 
ponds to about 1 mW of power. Compared with 
the power dissipation of the sensor chips, this is 
negligible and easily taken care of by the helium 
cooling. 



1 Muons are generated with an energy spectrum model- 
ing the one observed in MEG. 
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9.2 Vertex distribution 



The distances between tracks on the target can 
be reconstructed already online and used by the 
event filter farm to reject frames containing only 
background, see section [14) The only physics pro- 
cess exhibiting three tracks from the same vertex 
is the radiative muon decay with an internal con- 
version at a rate of 3.5 ■ 1CP 5 . 

The simulated distribution of vertices (more 
precisely: intersections of simulated particles with 
the target) is shown in Figures 9.2 and 9.3 In 



F 2500 



the longitudinal direction, the effect of the thicker 
material in the back part can be clearly discerned, 
whereas for the transverse view, the "shadows" of 
the target suspension are visible in the projected 
beam profile. 



Figure 9.4 shows the shortest distance between 
two vertices in a frame for 2 • 10 8 Hz muon stop 
rate. Less than 10% of the frames have tracks 
that come within 1 mm of each other on the tar- 
get. Figure [975] shows the shortest distance within 
which three tracks approach on the target sur- 
face; one of the tracks has to be assigned negative 
charge, either because it is a true electron or a re- 
curling positron track. Figures |9.6| show the same 
distributions for 2 ■ 10 9 Hz muon stop rate. Here 
all frames have two tracks approaching to closer 
than a millimeter, but a three track coincidence 
requirement still has a considerable suppression 
power, which can surpass a factor of 10~ 3 , if re- 
curlers can be identified with high efficiency (see 
also chapter 14 1. 



9.3 Alternative Designs 



9.3.1 Material Alternatives 

If the aluminium foil design proves unworkable or 
not mechanically stable enough, it could be re- 
placed by an equivalent design in carbon fibre re- 
inforced plastic (CFRP), where the material thick- 
ness would be approximately double. 200 um thin 
CFRP structures were built e.g. for the CMS pixel 



detector upgrade 125 



Another material option would be the use of a 
low density foam-like material such as Rohacell, 
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Figure 9.2: Vertex distribution along the beam 
direction. 
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Figure 9.3: Vertex distribution transverse to the 
beam direction. 



as used in the SINDRUM experiment 126 . 
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Figure 9.4: Shortest vertex- vertex distance inside 
a readout frame with 10 tracks on average (phase 
IB). Note that every crossing of a simulated elec- 
tron/positron track is counted as a vertex. The 
red histogram shows the contribution from in- 
ternal conversion decays. 
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Figure 9.6: Shortest vertex-vertex distance inside 
a readout frame with 100 tracks on average (phase 
II) . Note that every crossing of a simulated elec- 
tron/positron track is counted as a vertex. The 
red histogram shows the contribution from in- 
ternal conversion decays. 
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Figure 9.5: Shortest distance containing three Figure 9.7: Shortest distance containing three 
s + e _ e + inside a readout vertices consistent with e + e~e + inside a readout 

frame with 100 tracks on average (phase II). 
Note that every crossing of a simulated elec- 
tron/positron track is counted as a vertex; charge 
assignments are made purely on the apparent 
curvature, i.e. recurling positrons are counted as 
electrons. The red histogram shows the contribu- 
tion from internal conversion decays. 



frame with 10 tracks on average (phase IB). 
Note that every crossing of a simulated elec- 
tron/positron track is counted as a vertex; charge 
assignments are made purely on the apparent 
curvature, i.e. recurling positrons are counted as 
electrons. The red histogram shows the contribu- 
tion from internal conversion decays. 
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9.3.2 Active target 

We have also considered the use of our Kapton- 
sensor assemblies (see chapter [To]) as an active tar- 
get. This would lead to a vertex separation ability 
in the order of the pixel size (80 um) , as opposed 
to the w 200 um expected from track extrapola- 
tion. At an appropriate inclination angle, a wedge 
of two chips with 6x2 cm 2 active area and 50 um 
thickness presents a comparable amount of mater- 
ial to the beam (however, the cabling and cooling 
required would add significant additional material 



at least in the downstream direction). Such an 
arrangement would sacrifice ^-symmetry. 

The hit rate expected in an active target exceeds 
3 GHz for high intensity running, corresponding 
to about 8 KHz per pixel, or 500 MHz per reticle, 
far beyond the 80 MHz expected in the innermost 
sensor layer. This in turn would necessitate the 
development of a completely new read-out block 
for the sensor chips. The gain in selectivity when 
going from sa 200 um vertex resolution to 80 um 
does not justify the cost and technical risks asso- 
ciated with the active target. 
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Chapter 10 




The Mu3e Pixel Detector 



The Mu3e pixel tracker is to be built from High- 
Voltage Monolithic Active Pixel Sensors (HV- 
MAPS) thinned to 50 um. Signal and power 
lines are aluminum traces on a Kapton flex-print, 
which, together with a Kapton prism, also serves 
as a support structure. The detector should be 
cooled with gaseous helium. 



10.1 HV-Maps Sensor 

We propose to use Monolithic Active Pixel Sensors 
(MAPS) as tracking detectors as they integ- 
rate sensor and readout functionalities in the 
same device and thus greatly reduce the mater- 
ial budget. Classical concepts like hybrid designs 
have usually a higher material budget due to ad- 
ditional interconnects (bonds) and extra readout 
chips, which compromise the track reconstruction 
performance especially at low track momentum. 




N-well 




E field > 


P-substrate 






Figure 10.1: Sketch of the MAPS detector design 



from 127 



In the first MAPS designs ionization charges 
were collected mainly by diffusion with a tim- 
ing constant of several hundreds of nanoseconds. 



MAPS designs with high bias voltages exceeding 
50 V, however overcome this problem by collecting 
charges via drift and provide timing resolutions of 
better than 10 ns. 

We propose to use the High Voltage MAPS 
design with the amplifier electronics completely 
implemented inside the deep pixel N-well, which 



128 129 



was first proposed in 127 and since successfully 
tested 



see also section 10.4 
Figure [TO . 1 1 shows a sketch of a Monolithic Pixel 
Detector. The readout circuitry allows an efficient 
zero suppression of pixel information and the im- 
plementation of timcstamps to facilitate the as- 
signments of hits between different pixel layers. 

For readout designs providing 50 — 100 ns tim- 
ing resolutions power consumptions of about 



150mW/cm 2 are expected 130 . 

Because of the small size of the active deple- 
tion zone, the detectors can be thinned down to 
50 um or less. By thinning, the material budget 
can be significantly reduced and becomes, aver- 
aged over the tracking volume comparable to or- 
dinary gaseous detectors. 

A further advantage is that HV-MAPS can be 
implemented in a "cheap" commercial process. 
We use the AMS/IBM 180 nm HV-CMOS pro- 
cess 131 , which was developed mainly for the 
automotive industry, and thus offers long-term 
availability as well as being specified for a very 
wide range of operating conditions. The process 
offers a maximum reticle size of 2 x 2 cm 2 . One 
of the few disadvantages of the process is the fact 
that the first metalization layer is in copper, thus 
introducing a small amount of medium Z mater- 
ial. There are however plans for replacing also 
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Small Sensor Large 


Sensor 


Pixel Size [jim 2 ] 


80 x 80 I 


50 x 80 


Sensor Size [cm 2 ] 


1.1 x 2 


2x2 


Assembly 


1 x 3 


1 x 3 


Assembly size [cm 2 ] 


1.1 x 6 


2x6 


Max. LVDS links 


4 


2 


Bandwidth [Gbit/s] 


3.2 


1.6 


Table 10.1: 


Sensor specifications 





serialized on a 800Mbit/s low voltage digital sig- 
naling (LVDS) link. The sensor is configured via a 



that layer with an aluminium metalization in a 
future version of the process. 

Radiation-tolerance studies of the HV-MAPS 
sensors in 180 nm technology are ongoing also on 
other projects (e.g. ATLAS pixel R&D). Several 
test chips with similar pixel electronics as MUPIX 
have been irradiated at PS (CERN) up to doses 
between 80 and 430 MRad (the latest corresponds 
to a fluence of nearly 1 • 10 16 neq/cm 2 ). The res- 
ults are promising. Despite of the use of standard 
NMOS layouts, the chip irradiated to 80 MRad 
still detects the particles radiated by a 90 Sr source. 
The setup irradiated to 430 MRad is strongly ac- 
tivated, and no accurate detection of 90 Sr sig- 
nals is possible. The chip was able to detect 
particles from the beam up to 410 MRad dose. 
The main radiation effect is that the electronics 
suffers from ionizing effects such that it is difficult 
to find a proper operation point for the pixel amp- 
lifier. As the Mu3e experiment is performed at 
a muon beam-line the requirements on radiation 
hardness are not comparable to those at hadron 
colliders like the LHC. The radiation tests of the 
HV-MAPS sensors done so far all indicated that 
there will be no radiation damage even at highest 
muon rates at phase II. 



10.2 Sensor specification 

We plan to use two types of sensors in the Mu3e 
experiment, a smaller one for the inner layers and 
a larger one for the outer layers, see Table |10.1| 
The pixel size is 80 x 80 pm 2 , much smaller than 
the multiple scattering contribution. 

The wafers are to be thinned to 50 pm. If the 
yields permit it, we will cut strips of three sub- 
sequent sensors from the wafers and mount them 
in one piece. 

The sensor output is zero suppressed and con- 
sist of time-stamps and addresses of hit pixels, 



J TAG interface [132] , Including supply voltages, 
we expect about 30 pads (and thus bond wires) to 
connect the chip to the Kapton flex- print. 

10.3 Path towards the Full Sensor 

10.3.1 The MUPIX Prototypes 

First purpose-built sensor prototypes (the MUPIX 
series of chips) became available in 2011. 

MUPIX1 and 2 

The MUPIX 1 and 2 are small demonstration pro- 
totypes with a matrix of 42 x 36 pixels of 30 x 
39 pm 2 size for an active area of approximately 
1.8 mm 2 , see Figures 10.2 and 10.4 Each pixel 



consists of the sensor diode, a charge-sensitive 
amplifier and a source follower to drive the signal 
to the chip periphery. In addition there is a capa- 
city allowing to inject test charges. On the peri- 
phery, a comparator turns the analog signal into 
a digital time-over-threshold (ToT) signal. The 
threshold of the comparators is set globally for 
the chip and adjusted pixel-per-pixel with a 4-bit 
tune digital to analog converter (DAC). See Fig- 
ure 1 10 ,3| for an overview of the pixel electronics. 
In the test chips, the comparator output of an 
individual pixel can be observed via a dedicated 
output line. Alternatively, the whole chip can be 
read out via a shift register, where all the avail- 
able information is whether a particular pixel saw 
a signal during an active gate. 

The MUPIX 1 chip had an issue with feedback 
in the comparator that occasionally led to double 
pulses. This issue has been fixed in MUPIX 2, 
which in addition contains temperature sensors. 

MUPIX3 

In August 2012, we submitted the MUPIX 3 chip, 
a major step towards the final sensor. The new 
chip has 40 x 32 pixels of 92 x 80 pm 2 size for 
an active area of approximately 9.4 mm 2 , see Fig- 
ure |10.5| It implements the full digital column 
logic, allowing for address generation and serial 
readout of zero-suppressed data. In addition, 
MUPIX3 has faster signal shaping. 

The main differences with the final sensor are 
the lack of a high-speed LVDS output, buffers in 
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J 



Figure 10.2: Design view of the MUPIX2 chip (ac- 
tual size about 1.8 by 2.5 mm). 




Figure 10.3: Schematic of the pixel cell analog 
electronics in the MUPIX chips, nw stands for 
n-well. 
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Figure 10.4: MUPIX2 sensor on ceramic carrier 
with one Euro-cent for scale comparison. 



Si ili aiiu-; Si. SiiSV.Si u Si Si Si &) i.Si: iS'j3 

■ ■ 

' ' - 

.! . Si Si. -J ..3 Si : St .Si :Si .Si. , 

■ £ | : ( ^ i ■ ^ q o 

i n « iri =r n . '-(: .St iri n 




■-. - 

b ill I'l 




Figure 10.5: Design view of the MUPIX3 chip (ac- 
tual size about 4 by 5 mm) . 



the columns and the chip- wide hit collection logic. 
For this prototype, the corresponding logic will be 
emulated off-chip in an FPGA. Also several con- 
trol voltages which should be generated on-chip in 
the final version are currently produced on a test 
printed circuit board (PCB) in order to allow for 
easier debugging. 

We have just received the first MUPIX 3 
samples. First test results will become available 
early 2013. 



10.3.2 Plans for 2013 

As soon as the first results from the MUPIX 3 
sensor are available, we will prepare another 
multi-project-wafer (MPW) run, implementing 
the remaining digital logic and addressing poten- 
tial issues discovered with MUPIX 3. This should 
clear the path for an engineering run in the second 
half of 2013, opening up the possibility to build a 
full scale tracker prototype. 
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10.4 Characterization of the Proto- 
types 

Wc have studied the properties of the MUPIX 1 
and 2 prototypes using injection pulses, LEDs, 
laser diodes, X-rays, radioactive sources and test 
beam measurements. In the following, we will out- 
line the core results; details of the findings can be 
found in a master 133 and a bachelor thesis |134] . 

The test setup was based on a test board hous- 
ing the chip itself and a logic interface card. The 
chip test board has voltage regulators for the sup- 
ply voltages, digital to analog converters for the 
threshold and injection pulse height and a flat rib- 
bon connector to the logic interface card. This in- 
terface card uses a FPGA both to program and 
readout the MUPIX chip and to communicate 
with a PC via USB. These two boards together 
with a C++ control software have been used for 
all tests described here. 

10.4.1 Signal to Noise Ratio 

Injection pulses can simulate the charge depos- 
ition of ionizing particles. A scan of the injection 
pulse height at constant threshold would deliver 
a step function if no noise were present. Due to 
the additional noise on top of the injection pulse 
the step function is transformed to an error func- 
tion with a finite slope. In turn the width of the 
step function can be used to determine the noise 
as a function of the applied threshold. A 55 Fe 
source has been used as reference signal for the 
determination of the signal to noise ratio (SNR). 
The injection pulse height corresponding to the 
55 Fe signal is found by matching the time over 
threshold of both signals at a given threshold. For 
the MUPIX 2 chip the signal to noise ratio as 
a function of threshold ranges between 21.5 and 
35.7, see Figure [T0.6| Any signal to noise ratio >9 
is considered good. 

10.4.2 Pixel to Pixel Uniformity 

Using a similar method as for the measurement of 
the signal to noise ratio, the pixel to pixel uni- 
formity can be determined and optimized. At 
a constant injection pulse height a scan of the 
threshold setting for the entire pixel matrix is car- 
ried out. For each pixel the error function fit is 
performed and the 50 % value of the error func- 




0.9 0.92 

Threshold [V] 



Figure 10.6: Signal to noise rate for the MUPIX 2 
chip. The signal size is taken from measurements 
with a 55 Fe source, whilst the noise is measured 



in a injection pulse scan, sec 133 for details. 



be seen in Figure 10. 7( ;) the distribution of the 
threshold values for the pixels of the same chip 
is quite broad. In order to compensate for these 
inequalities, a threshold offset, the so called tune 
DAC (TDAC) value can be set for each pixel. The 
best TDAC values for the pixels are found by run- 
ning an automated procedure, see 133 . Figure 



tion rising edge plotted, see Figure 10.7 1). As can 



10.7 d) and d) show that the threshold spread (re- 
flecting the response uniformity) after tuning the 
threshold offsets has decreased by almost an order 
of magnitude from 19.6 mV to 2.3 mV RMS. 

10.4.3 Pixel Response Time 

The timing characteristics of the MUPIX chip is 
very important because of the requirement to run 
at high muon decay rates. A fast response of the 
MUPIX chip leads to precise timing of hits and 
helps suppressing combinatorial background. In 
the inner detector layers the double pulse resolu- 
tion should be good as the hit rate is up to 2 kHz 
for each pixel. 

Timing studies have been carried out using a 
LED driven by a pulse generator to stimulate the 
sensor. The discriminator output of a single pixel 
is then compared to the second output of the 
pulse generator on an oscilloscope. In this con- 
figuration it is possible to measure the latency 
between generator and pixel output and the time 
over threshold (ToT) of the pixel. By plotting 
latency and ToT for many threshold values one 
obtains the pulse shape, see Figure |10.8| The 
measured pulse shape reflects the RC-CR shaping 
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Figure 10.7: Uniformity of the pixel response before and after tuning of the pixel DAC values, from 133 
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Figure 10.8: Measured pulse shape of a MUPIX 
pixel in response to an LED pulse, fitted with the 



expectation from RC — CR shaping, from 133 
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Figure 10.10: Time-over-threshold (corresponding 
to energy) spectrum of a 55 Fe radioactive source, 
from [133]. 
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Figure 10.9: Double pulse resolution of the 



MUPIX2 chip, from 133 



characteristic of the pixel electronics. The charge 
collection process in the depletion zone is much 
faster and has little influence on the measured tim- 
ing. 

By generating double pulses and dividing the 
number of unresolved pulses at the pixel output 
by the number of all pulses, the double pulse res- 



olution can be determined. Figure 10.9 shows this 
ratio as a function of pulse to pulse delay. An er- 
ror function has been fitted and the point of 1 % 
un-resolved pulses gives a double pulse resolution 
of 2.7 us consistent with the expectation from the 
chosen shaper characteristics. 



10.4.4 Measurements 
Sources 



with Radioactive 



Radioactive sources allow to test the MUPIX 
chips with real particles. For a 55 Fe source the 
time over threshold (ToT) distribution, which is 
an measure for the energy distribution, is shown 
in Figure |10.10| X-ray fluorescence in the range 
of 4 to 18keV has been used to derive a relative 
energy resolution of 10 to 20 %. 

The 55 Fe peak in the ToT spectrum has been 
used to study the influence high voltage (HV) and 
temperature on the pixel signal. Figure |1 . 11 1 
shows the 55 Fe peak as a function of the HV, 
revealing first a steep rise as the depletion zone 
grows. Between 20 and 60 V, the ToT stays al- 
most constant. At even higher voltages the elec- 
tric field in the depletion zone becomes strong 
enough for the creation of secondary electron- hole 
pairs, which results in a signal amplification. The 
temperature dependence of the pixel sensor is im- 
portant, because in the later experiment the differ- 
ent parts of the detector will operate at different 
temperatures varying by 30-40 °C. In Figure [T"0.12| 
the 55 Fe peak of the ToT spectrum is plotted as a 
function of temperature between 20 to 60 °C. As 
can be seen the ToT goes down from 1.6 to below 
1 us, with a slope of 16.3ns/°C. Measurements 
with test pulse injection and simulation studies 
confirm that the observed temperature depend- 
ence is a feature of the charge sensitive amplifier 
of the pixel analog stage. 



48 




An Experiment to Search for the Decay fi — > eee 



1.5 — 
1.4 — 



1.2 - 
L.I — 



70 80 90 

High Voltage [V] 



Figure 10.11: Position of the Fe peak in depend 



ence of the applied high voltage, from 133 
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10.4.5 Testbeam results 

In August 2012 the MUPIX2 chip was tested at 
the CERN SPS. The beam from the SPS was a 
170 GeV pion beam, chosen for little multiple scat- 
tering effects due to the device under test (DUT) . 
In two separate data taking periods the MUPIX2 
was mounted inside the TimePix telescope |135| . 
This silicon pixel telescope has four layers before 
and four layers after the DUT, providing very pre- 
cise pointing resolution of 5 pm. The MUPIX2 
chip was tested facing the beam and under an 
angle of 45°. For both beam periods a threshold 
scan was performed in order to derive the effi- 
ciency as a function of threshold, but no threshold 
calibration to equalize the gain of the pixels was 
performed. Figure |10.13| shows the resolution in 
x and y direction of the MUPIX2 chip facing the 
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Figure 10.13: Resolution of the MUPIX 2 chip 
measured in a testbeam with 170 GeV pions and 
the TIMEPIX telescope. The measured resolu- 
tions correspond very well to the expectation from 
the 5 um resolution of the telescope and the pixel 
pitch in x and y direction, respectively. 




Figure 10.14: Tool for layer 3 segment assembly. 



beam and tested at three different threshold val- 
ues. The measured resolution of 11.2 pm in x and 
15.4 pm in y corresponds to the expected resolu- 
tion given the telescope resolution and the pixel 
size. Further results are still subject to studies 
and will follow soon. 



10.5 Mechanics 

The pixel detector mechanics has been optimized 
for very low material budget in the active detector 
region. Additional requirements are mechanical 
stability, resistance to temperatures over a wide 
range and a modular design for ease of assembly 
and repair. It is proposed to build the frame for 
the silicon pixel detector from thin Kapton foil. 
The sensors are glued and bonded on a flex-print 
and then mounted onto the Kapton frame. 
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Figure 10.15: Tool for layer 3 segment assembly. 



Figure 10.18: Mounting tool for the central pixel 
detector. 




Figure 10.16: Mechanics of the central pixel de- 
tector 




Figure 10.17: Segmentation of the central outer 
layers. 



The Kapton foil used for the mechanical frame 
construction is 25 jim thin. It gains in mechan- 
ical stability as it is folded around a prism-shaped 
template and glued to plastic end-pieces. For this 
process tools have been built, which support the 
detector modules until they are ready for full sta- 
tion assembly, see Figures |10.14| and |10.15| In a 
separate tool the 50 um thin pixel chips are glued 
on the single layer Kapton flex-print. Alignemcnt 
groves secure the correct position of the chips. In 
the next step the chips are wire bonded to the 
flex print. With a vacuum lift tool the hex-prints 
are then positioned and glued to the mechanical 
frames. The radiation length of a pixel detector 
layer is summarized in Table |10.2| 

The pixel detector is build from layers of four 
different sizes and prism shapes, see Figure [T0.16| 
As the digital readout circuits of the pixel chips 
create an approximately 0.5 mm wide dead area, 
there is a 1 mm overlap to the adjacent sensor. 
The inner double layers have 12 cm active length. 
Layer 1 has 12 and layer 2 18 sides of 1 cm 
width. Each inner layer is assembled from two 
half-modules. As a consequence the plastic end- 
pieces are half moon shaped. When all four 
half modules of the inner detector layers are pro- 
duced and tested they are mounted to two thin 
rim wheels. A mechanical prototype for the in- 
ner double layer has been constructed from 25 um 
Kapton foil both for the frame and the flex-print 
layer, while the pixel chips have been simulated 
with 100 um thick glass plates. Glass of 100 um 
thickness is of comparable flexibility as thinned 
silicon, which is shown in Figure |10.19| The res- 
ulting mechanical unit is surprisingly sturdy and 
fully self supporting, as can be seen in Figure 
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Component 


THICKNESS 

[urn] 


X/X 

[%] 


Kapton frame 


25 


0.018 


Kapton flex-print 


25 


0.018 


Aluminum traces 


15/2 


0.008 


(50 % coverage) 






HV-MAPS 


50 


0.053 


Adhesive 


10 


0.003 


Full detector layer 


125 


0.100 



Table 10.2: The pixel detector layer radiation 
length is dominated by the HV-MAPS chips. 



|10.20| Layer 3 and layer 4 have a three times lar- 
ger active length of 36 cm. The sides of these outer 
layers are 19 mm wide. The layers 3 and 4 have 
24 and 28 sides. An outer double layer module 
combines four sides, so layer 3 consists of 6 and 
layer 4 consists of 7 modules, see Figure |10.17| 
The station assembly is done in a special mount- 
ing frame, combining the inner two layers with 
the modules of the outer two layers on large rim 



wheels (Figure 10.18). Re-curl stations will be as- 
sembled accordingly. 

10.6 Cooling 

The cooling system for the pixel detector must be 
capable of keeping the temperatures at a reason- 
able level (< 70 °C) and should at the same time 
add very little extra material to the active volume 
of the detector. Cooling with gaseous helium has 
been chosen as it appears to offer a reasonable 
compromise between cooling potential and radi- 
ation length. 

Given its size, the heat load introduced by the 
pixel detector is considerable. Table |10.3| de- 
tails the contributions of the different parts of 
the pixel detector assuming a realistic heat load 
of 150mW/cm 2 . It is worth mentioning that 
the heat load is strongly dependent on the clock 
speed at which the pixel chips are running. Pixel 
chips with very low occupancy, for example those 
installed in the recurl stations, could optionally 
be operated at half the clock speed which would 
largely reduce the total thermal power of the pixel 
detector. 

In order to keep the temperature of the coolant 
equally low in all parts of the detector, it is fore- 



Detector Part 


size [cm 2 ] 


power [W] 


layer 1 


158.4 


23 


layer 2 


237.6 


35 


vertex layers (1+2) 


396 


59 


layer 3 


1728 


259 


layer 4 


2016 


302 


outer layers (3+4) 


3744 


561 


central detector (1-4) 


4140 


621 


detector with 






2 recurl stations 


11628 


1744 


full detector 


19116 


2867 



Table 10.3: 
150mW/cm 2 



Pixel detector heat load for 



seen to place nozzles between the support struc- 
tures pointing at different parts of the sensor lay- 
ers. 

In order to study the potential of gaseous he- 
lium cooling a bachelor thesis has been carried 
In this thesis a thin sandwich of alu- 
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out 

minum foil, Kapton and silicon was heated induct- 
ively and cooled with a helium flow. For the 8.5 cm 
long sample heated with 100mW/cm 2 a flow of 
0.4 m/s was necessary to maintain a temperature 
difference between helium and sensor of 32 K. Fur- 
ther studies using a full scale central detector pro- 
totype in a wind tunnel are on the way. This full 
scale prototype will be assembled from aluminized 
Kapton foil such that Ohmic heating can be ap- 
plied. It is planned to run the upcoming tests with 
the helium heat exchanger acquired for phase I of 
the experiment. 



10.7 Alternative Technologies 

Since the Mu3e pixel detector plays the key role in 
the experiment, it is important to carefully com- 
pare the chosen HV-MAPS technology with al- 
ternatives. To do so one has to have the detector 
requirements in mind, see chapter [5j The goal of a 
branching ratio sensitivity of 10 -16 leads to a high 
muon decay rate on target of 2 • 10 9 Hz. In order 
to suppress combinatorial background a very good 
vertex resolution of O(200 pm) is needed. The 
suppression of the background from the conver- 
sion decay /i — > eeevv requires very precise mo- 
mentum reconstruction of better than 0.5 MeV, 
for the 10~ 16 branching ratio sensitivity. The 
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momentum resolution for the low energetic decay 
products of 10 MeV to 53 MeV will be governed by 
the multiple scattering in the detector. Only con- 
cepts with extremely low radiation length of < 1 % 
in total can achieve this. 

In the following different technology options 
for the MuSe tracking detector are compared. 
The result of this comparison is summarized in 
table HUH 

10.7.1 Time-Projection Chamber 

Time projection chambers are a very attractive 
concept for tracking detectors with good spacial 
resolution and very low radiation length. The ex- 
ample of the ALICE TPC |137| has demonstrated 
that very high track densities can be handled 
with such a detector. Over the last two dec- 
ades many studies for TPCs with very low radi- 
ation length, i.e. 0.12 % for a proposed chamber at 
SLAC |138] , have been carried out. The drift time 
of an electron in a possible TPC for the Mu3e ex- 
periment with 2 m total length would be in the 
order of 50 us, using the proposed gas mixture 
of Helium: C0 2 :Isobutane (83:10:7%). This time 
of 50 us at 2 • 10 9 Hz decay rate means that the 
DAQ and event filter farm would have to handle 
frames containing >10 5 muon decays, which is not 
realistic. Another problem would be strong space 
charge effects due to continuous very high track 
density, see |139| . These two aspects rule out the 
usage of a TPC in the phase IB and phase II of 
Mu3e, while it would be possible to use a TPC in 
phase IA. 

10.7.2 Other Pixel Technologies 

The usage of pixel detectors in particle and heavy- 
ion physics has recently lead to a variety of thin 
detectors capable of dealing with high particle 
fluxes. It is worth to discuss the solutions de- 
veloped for the upgrades of the STAR, BELLE 
and LHCb vertex detectors. 

The STAR upgrade vertex detector [140 
is based on Monolithic Active Pixel Sensors 
(MA PS), more precisely the MIMOSA 28 sensor 
|141| . The vertex detector under construction is 
located around the interaction point at r=2.5cm, 
8 cm and has an intrinsic hit resolution of <6 um. 
An ultralight carbon fibre support structure has 
been developed with a thickness of 200 um. The 
radiation length is X/Xo=0.079 % per layer in the 
active region. The power dissipation is about 



170mW/cm 2 and forced air cooling at lOm/s is 
used [142; . The charge collection for the MIMOSA 
chips is based on electron diffusion, resulting in 
charge collection times of 50 ns. This has to be 
compared to the charge collection times of drift 
based silicon detectors like silicon strip and hy- 
brid silicon pixel detectors of a few ns. The signal 
integration of the MIMOSA 28 chip is 185.6 us. 
This long integration time corresponding to >10 5 
decay events rules out the usage of the MAPS de- 
veloped for STAR vertex upgrade in Mu3e, 

Bellc2 is building a lightweight silicon vertex 
tracker based on DEPFET pixel sensors 43 . The 



Belle2 silicon vertex detector has two layers at 
r=1.4cm, 2.2 cm. The monolithic sensors are 
only 75 um thick and self supporting, resulting 
in a radiation length of X/Xo=0.18% per layer 
143 . The cooling system uses liquid CO2 for the 



readout electronics and cold dry air with forced 
convection for the sensors. The sensor chip is 
equipped with an integrated read out amplifica- 
tion, the charge is then accumulated internally. 
The readout is enabled externally and done in a 
row wise rolling shutter mode. The readout time 
for the entire chip is 20 us. As for the MIMOSA 28 
sensor the relatively long readout time makes the 
Belle2 DEPFET sensor unsuitable for the Mu3e 
detector. On top of this the Belle2 DEPFET 
sensor produces analog output signals and re- 
quires further electronics for analog to digital con- 
version, baseline subtraction, discrimination and 
zero-suppression . 

The proposed pixel detector for the LHCb ver- 
tex upgrade belongs to the TIMEPIX family of 
chips. The main goal of the LHCb upgrade is to 
move from a 1 MHz readout to a 40 MHz readout 
in order to make a track based LO-trigger possible. 
The power consumption per chip is 1.5W/cm 2 



144 , which is quite high. The sensor thickness is 



150 um to 200 um. The LHCb Velo upgrade pixel 
sensor is a hybrid design, so in addition to the 
sensor there is also the readout chip thickness, the 
solder bump bonds, aluminum traces and cooling 
elements summing up to X/X o =0.73% per layer. 
In LHCb the vertex detector is inside a vacuum 
tank around the interaction point, so the cool- 
ing scheme relies on a high thermally conductive 
spine (diamond) in combination with a liquid CO2 
cooling outside the active area. Alternative cool- 
ing studies look at micro-channel cooling which 
would extend the CO2 cooling to the pixel chips 
in the active area. The pixel sensor is designed to 
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Parameter 


required 


Mu3e 
HV-MAPS 


He TPC 


STAR 
MIMOSA28 


BELLE 
DEPFET 


LHCb 
VeloPix 


sensor thickness X/Xo [%] 
readout cycle [us] 
power [mW/cm 2 ] 


<0.1 

< 0.1 

< 200 


0.075 
0.05 
150 


50 


0.079 
185.6 
170 


0.18 
20 
100 


0.73 
0.025 
1500 



Table 10.4: Alternative technology comparison, only the crucial parameters are listed. 



have less than 25 ns time-walk at 1000 electrons. 
The output bandwidth is > 12 Gbit/s with on-chip 
zero suppression. The readout is hit driven and 
asynchronous. While the readout speed would be 
ideal for the Mu3e detector, the radiation length 
per layer of X/Xo=0.73 % is too high. 

10.7.3 Conclusion 

In the field of tracking detectors with good res- 
olution, low radiation length and high rate cap- 
abilities there has been a very strong progress in 
the recent years, see Table |10.4| Time projection 
chambers with helium gas mixtures would provide 
excellent momentum resolution. Because of the 
very long drift times in a large TPC and the con- 
stantly high muon decay rates, it is not possible to 
do online track reconstruction in the event filter 



farm for decay rates exceeding 1 • 10 7 Hz, which 
corresponds to phase I A of the Mu3e running. As 
a change of technology between phases IA and IB 
is not desirable, the TPC option can only be con- 
sidered as a back-up solution. 

The above discussed pixel alternatives can be 
grouped in two categories. The systems based 
on the MIMOSA chip for the STAR upgrade and 
the DEPFET chip for Belle 2 have little radiation 
length of X/X o =0.079 % and 0.18 % per layer, but 
have long readout cycles. On the other hand the 
VeloPix based tracker for the LHCb upgrade can 
be read out at 40 MHz, but has a relatively large 
radiation length of X/Xo=0.73% per layer. In 
all three cases further R&D would be required to 
transform the existing designs into one fulfilling 
the requirements of all phases of the Mu3e exper- 
iment. 
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Chapter 11 




The Mu3e Fibre Detector 



11.1 The time of flight detector 




Figure 11.1: Schematic side view (left) and cross 
section (right) of the detector components. The 
Sci-Fi system is highlighted in blue. 

A cylindrical time of flight (ToF) detector com- 
plements the central silicon tracking system. It 
consists of a scintillating fibre (Sci-Fi) hodoscope 
with a radius of 6 cm and a length of 36 cm. The 
expected time resolution is of several 100 ps and 
a detection efficiency close to 100%. The main 
purpose of the ToF system is to measure very pre- 
cisely the arrival time of particles in order to allow 
for the matching with hits detected in the silicon 
detectors. This will help to reject pile-up events 
(accidental backgrounds) and allow for a charge 
(direction of propagation) measurement for recur 1- 
ing tracks. The ToF system will operate at very 
high particle rates up to several MHz per channel. 

A detailed R&D program is ongoing to prove 
the feasibility of the Sci-Fi detector and to help 
optimizing the design of the Sci-Fi detector. The 
R&D activity covers all the aspects of the ToF 



detector development: scintillating fibres, SiPMs, 
amplifiers, and readout electronics. The chal- 
lenging aspects of the ToF system are (i) the 
high rates per Sci-Fi readout channel and (ii) 
the large data flow generated by the readout di- 
gitizing electronics, which has to be handled in 
real time. Time resolutions of about 200 to 
300 ps have been already achieved with Sci-Fi 
hodoscopes with single-ended readout using multi- 
anode PMTs 145 . We aim at a similar time res- 



olution. 

In the current design, the scintillating fibre de- 
tector is composed of 24 Sci-Fi ribbons, each 36 cm 
long and 16 mm wide, as illustrated in Figure [Tl. 2 1 
The ribbons will be supported by a carbon fibre 
mechanical structure. The use of 250 um dia- 
meter multi-cladding scintillating fibres produced 
by Kuraray is envisaged. Three to five Sci-Fi 



layers (Figure 11.2 right) are staggered such as 
to minimize empty spaces between fibres. Fig- 
ure 111.31 shows a Sci-Fi ribbon and its cross sec- 
tion. In the figure the uniform staggering of the 
Sci-Fi layers can be seen. The overall thickness of 
the Sci-Fi arrays varies depending on the number 
of layers between 0.6 to 1 mm. This thickness has 
to be kept as low as possible, to reduce the mul- 
tiple scattering to a minimum compatible with the 
required performance (i.e. time resolution and ef- 
ficiency). The empty space between adjacent rib- 
bons will be minimized in order to guarantee con- 
tinuous coverage. 

The scintillating light produced in the fibres will 
be detected with SiPM arrays at both fibre ends. 
The choice of SiPMs as a photo-detection device is 
based on the fact that they are very compact ob- 
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Figure 11.2: Left: Overall view of the Sci-Fi barrel ToF detector. The diameter of the detector is about 
12 cm, the length about 36 cm. Middle: Detail showing the SiPM ceramic supports and the Sci-Fi arrays. 
Right: Details of the Sci-Fi ribbons. In this figure, the ribbons consist of five staggered Sci-Fi layers. 




Figure 11.3: Photos of a Sci-Fi ribbon (left) and of its cross section (right). The uniform staggering of 
4 Sci-Fi layers is clearly visible. 



jects that can be operated in high magnetic fields 
with high gain (~ 10 6 ) and at high counting rates. 

Two different readout schemes for the fibres are 
being considered and investigated: 1) fibres are 
grouped in vertical columns (this is implemented 
by the structure of the readout photo-device) , 2) 
each fibre is read out individually. In the first 
design we envisage to use the 32 channel SiPM ar- 
rays with 50 x 50 jim 2 pixels available from Hama- 
matsu 146 147 . The active size of the sensor is 
8 x 1 . 1 mm 2 with a total surface of 9 mm 2 . The 
pixels are arranged in columns, corresponding to 
an effective pitch of 250 um. Two such SiPM 
sensors will be assembled side by side with al- 
most no dead region, giving a 16 mm wide photo- 
sensitive region with 64 readout channels, and 
matching precisely the width of the Sci-Fi ribbons. 
The photo-detectors will be directly coupled to the 
Sci-Fi arrays to maximize the light collection effi- 
ciency, and consequently the timing performance. 
In total 96 SiPM arrays with 32 channels each will 
be required to readout the Sci-Fi tracker at each 



end, for a total of 2 x 24 x 64 (~ 3000) readout 
channels. 

A detailed understanding of the expected occu- 
pancies in the the Sci-Fi detector is required. The 
rates depend also on the background in the de- 
tector generated by secondary interactions in the 
materials of the different components. A particle 
crossing the Sci-Fi ribbons at 90°, will excite one 
or two adjacent Sci-Fi columns, as sketched in Fig- 
11. 4| A particle crossing the same Sci-Fi rib- 



bon, for example at 45°, will excite 4 to 5 adja- 
cent Sci-Fi columns for a 5 layer ribbon and 2 to 
3 for a 3 layer ribbon, thus increasing significantly 
the overall occupancy of the detector. At an ex- 
pected rate of 1.5 • 10 9 muon decays per second 
using 1500 Sci-Fi channels, the estimated rate in 
the Sci-Fi system is 5 MHz on average. At the 
same time the signal will be spread out over sev- 
eral SiPM channels. A well designed clustering 
algorithm will be required to reconstruct the tim- 
ing and position of the crossing particle. 

The lowest possible detector readout occupancy 
can be achieved by reading out each fibre indi- 
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vidually. In absence of dark current noise and 
background, the occupancy is given by the rate 
of detected electrons from muon decays and is ex- 
pected to be 1MHz per fibre. When reading out 
each multi-cladding fibre individually, one has to 
take into account that the scintillating light travels 
preferentially in the cladding of the fibre and exits 
the fibre in a cone with an aperture of about 45°. 
In order to collect all the scintillating light exit- 
ing the fibre, and thus ensure high detection effi- 
ciency, the photo-detector has to be wider then the 
fibre diameter, at least 100 pm around the fibre. 
This compensates also the misalignment, of the or- 
der of few tens of microns, between the fibres and 
the SiPMs. In order to avoid the light from one 
fibre to spill over to a neighboring photo-detector 
channel (optical cross-talk) the fibres and photo- 
detectors have to be separated by 500 pm center- 
to-center. To summarize, the estimated, optimal 
size of the photo-detector for reading out each 
250 pm fibre individually is 400 x 400 pm 2 (for a 
square device), and a center-to-center separation 
of 500 pm. 

A transition region between the highly compact 
Sci-Fi ribbons and the photo-detectors of about 
5 cm will be required to fan-out the Sci-Fis to the 
desired center-to-center separation of 500 pm be- 
fore coupling the fibres to the SiPMs. To assure 
the precise alignment, the fibres will be glued in 
sockets with a precision of 10 pm. The SiPMs will 
be coupled to the Sci-Fis using the same sockets. 
Reading out each fibre individually will require a 
total of about 9000 readout channels. 

In case the SiPMs could not be coupled dir- 
ectly to the Sci-Fis (because of mechanical con- 
straints inside the limited space), optical fibres 
coupled to the Sci-Fis will be used to transmit the 
scintillating light outside of the solenoid, where 
they will be coupled to the same type of SiPM ar- 
rays, however with some loss of light. In this case 
the use of multi-anode PMTs or multi-channel 
plate PMTs could be also envisaged as alternat- 
ive photo-detector device to the SiPMs for the 
readout of the Sci-Fi ribbons. 

Different manufacturers are also being investig- 
ated. The reported photo-detection efficiency of 
the Hamamatsu SiPMs is not that high (PDE ~ 
25 %) . Since the time resolution depends on the 
number of detected scintillating photons (<t t ~ 
l/yYiph) in order to obtain the best timing per- 
formances we need to maximize the PDE of the 
SiPM arrays. Discussions with Hamamatsu to ex- 




Figure 11.4: Sci-Fi occupancy vs. particle crossing 
angle and position. Dark fibres indicate higher 
energy deposit. 

plore various possibilities to increase the PDE are 
ongoing. We will also discuss with other compan- 
ies capable of producing such devices (we already 
contacted FBK in Italy and Ketek in Germany). 
Very likely it will be required to design SiPM 
sensors matching precisely the requirements of the 
Sci-Fi detector rather than using commercially 
available SiPMs. 

Performance studies of Sci-Fi ribbons coupled 
to SiPM arrays will be carried out in test beams 
at PSI. The test beam activities will include the 
study of the time resolution, rate capabilities, de- 
tection efficiency, tracking resolution, and uni- 
formity of the detector. Different readout elec- 
tronics will be also tested. 

11.2 Readout of photon detectors 

For the readout of the SiPMs different options are 
being investigated. Details are given in chapter 

EH 

One possibility is to use the well-established 
waveform digitizing technology based on the DRS 
switched capacitor array chip developed at PSI. 
The advantage of this technology compared to 
traditional constant fraction discriminators and 
TDCs is that pile-up can be effectively recog- 
nized and corrected for. Due to the bulkiness of 
this readout electronics, the digitizing electronics 
cannot be located in the proximity of the Sci-Fi 
detector. To transmit the electrical pulses from 
the SiPMs to the digitizing electronics two op- 
tions are being considered. In a first version the 
signals are transmitted without amplification on 
shielded low attenuation coaxial cables. Given 
the low amplitude of the pulses ( 1.5-2 mV per 
photo-electron) and potentially high noise from 
the digital readout electronics, the amplification 
of SiPM signals might be required. In this second 
version the pulses are first amplified about a factor 
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Figure 11.5: Simulated photon output integrated 
over 1000 positron crossing events for a specific 
fibre geometry (top) with the corresponding time 
distribution of the SiPM response using a constant 
fraction discriminator algorithm with a FWHM of 
about 300 ps (bottom). 



of ten and then transmitted. Hybrid amplifier 
circuits will be installed on the other side of the 
SiPM mountings. We already developed fast amp- 
lifiers using transistors (rise time ~ Ins, decay 
time ~ 10 ns) matching the input characteristics 
of SiPMs for optimal time performance and high 
rates. 

A second possibility is to use the STiC chip, 
which includes a fast discriminator and a TDC 
digitizing the time information. The advantage of 
this solution is the compactness of the chip that 
could be installed very close to the SiPMs, with 
no need to amplify and to transmit analog signals 



outside of the solenoid. A third option could be 
offered by time to digital converters implemented 
in FPGAs. 

11.3 GEANT simulations 

To optimize the overall design of the detector and 
of the Sci-Fi sub-system in particular, extensive 
simulation studies are being carried out. We are 
exploring the propagation of scintillating light in 
the fibres in a dedicated simulation that allows 
to compare the optical properties for any possible 
fibre and ribbon geometry. Meanwhile a complete 
set of simulation tools is available which allows a 
profound analysis of the different possible layouts 
of the Sci-Fi system. 

Ongoing simulation studies will help us to op- 
timize and work out the details of the final Sci-Fi 
ToF design. Given that a crucial requirement of 
the experiment is to minimize all involved mater- 
ials also a layer of reflecting paint or glue will add 
additional non active material to the detector and 
thus increase the overall material budget seen by 
the particles crossing the Sci-Fi ribbons, resulting 
in a deterioration of the overall momentum res- 
olution of the tracking system. This and similar 
effects can now be studied easily with the available 
tools. 

In combination with a simulation of the SiPM 
response |148| , we are now able to estimate 
whether the aimed time resolution of a few hun- 
dred picoseconds is feasible with the fibre geo- 
metry and the SiPMs under investigation. First 
studies show that the resolution can be achieved 
with the baseline design (Figure 11.5 1. Further 



fibre geometries and readout combinations are 
currently being evaluated. In addition to the sim- 
ulation test setups are being developed consist- 
ing of a complete ribbon with the SiPM detect- 
ors mounted on both ends to verify the simulation 
results and to get the prove whether the desired 
geometry can be constructed with reasonable ef- 
fort. Furthermore the test stand will be used to 
evaluate readout electronics. 

An additional source of concern is the cross talk 
and after pulsing observed in most APDs operated 
in Geiger mode. Both effects will lead to an in- 
crease of the occupancy of our electronics channels 
and thus have to be minimized. Since the Sci-Fi 
ribbons are read-out on both ends, this uncorrel- 
ated source of noise, in principle, can be rejected. 
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Chapter 12 




The Mu3e Tile Detector 




Figure 12.1: Drawing of a tile detector module 
with 48 x 48 tiles. 



The second component of the Time-of-Flight 
hodoscope is a detector consisting of scintillating 
tiles which are, as in the case of the Sci-Fi tracker, 
read out by SiPMs. The tile detector is located in 
the recurl station on the inside of the recurl pixel 



layer (see Figure 6.5 1. The detector aims at a time 



resolution of below 100 ps and an efficiency close 
to 100% in order to effectively identify a coincid- 
ent signal from three electrons and suppress com- 
binatoric background. The main challenge of the 
detector design is to achieve these requirements 
under the high hit rate which is expected, espe- 
cially in phase II. 

12.1 Detector Design 

The detector is divided into four modules (one for 
each recurl station) with a length of 36 cm and 



diameter of sa 12 cm. Each module is segmented 
into small tiles with a size of 0(1 cm 3 ), which are 
individually wrapped in reflective foil in order to 
increase the light yield. 

In the current design, each module consists of 
48 rings with a thickness of 5 mm, and each ring 



consists of 48 tiles (see Figure 12.1 1. This results 



in a total number of 2304 tiles per station, and 
a tile geometry of roughly 7.5 x 7.5 x 5 mm. For 
the outer two recurl stations, the number of tiles 
might be decreased to 36 x 36 due to the lower hit 



rate (see Figure 12.2 ), in order to reduce the costs. 
This detector geometry is the result of extensive 
simulation studies which will be continued in order 
to further optimize the detector performance. 

The tiles are made of plastic scintillator mater- 
ial which provides a fast light response. The main 
requirements for the scintillator are a high light 
yield, fast rise and decay time and a scintillation 
spectrum which approximately matches the spec- 
tral sensitivity of the SiPM. Different scintillators 
have been compared using a Geant4 simulation 
in combination with a simulation of the SiPM re- 
sponse. The best suited scintillator was found to 
be BC420 from Bicron. However, the differences 
between the individual scintillating materials were 
found to be small, and thus also more inexpensive 
alternatives are feasible. 

The scintillation light of each tile is read out by 
a SiPM which is directly attached on the inside of 
the tile. SiPMs are well suited for this application 
due to the compact size, insensitivity to magnetic 
fields, high photon detection efficiency (PDE) and 
excellent timing properties. The most important 
SiPM characteristics for this application are a high 
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The optical properties of a tile are parameterized 
using a separate Geant4 simulation of a single tile. 
The SiPM simulation GosSiP 11481 was used to 
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Figure 12.2: Hit rate per tile for phase II. The rate 
for phase I is a factor 10 to 20 smaller. 

PDE and gain in order to obtain a large signal 
and therefore a good time resolution. Further- 
more, a small signal decay time is desired in order 
to reduce signal pile up and thus increasing the 
signal efficiency. The S10362-33-050 MPPC from 
Hamamatsu with an active area of 3 x 3 mm 2 and 
3600 pixels is a suitable sensor for this applica- 
tion. A possible alternative sensor is the PM3350 
from Ketek which offers a larger PDE and gain 
compared to the Hamamatsu device, but has a 
longer decay time. A further option is the SensL 
MicroFB series, which will be available in the first 
quarter of 2013. These sensors are a very prom- 
ising alternative, as the devices will have an addi- 
tional fast output providing signals with a width 
of O(ns) and provide a high PDE and gain. Meas- 
urements and simulation studies comparing the 
different sensor types will be continued, in order 
to find the devices which is best suited for this 
application. 

For the readout of the SiPMs, two different op- 
tions are currently discussed. One solution is to 
digitize the SiPM waveforms using the DRS chip. 
The waveforms then have to be further processed 
externally, in order to determine the time-stamps. 
The second option is the STiC chip which includes 
a fast discriminator and a TDC digitizing the time 
information. Details on both chips can be found 
in chapter [13] 

12.2 Simulation 

Simulation studies using the full detector simula- 
tion described in section H~5l have been carried out 
in order to determine the detector performance. 



generate the signal waveforms from the hit inform- 
ation of the full detector simulation. The simu- 
lation has been done for a detector design with 
48 x 48 5 mm thick tiles per recurl station with 
BC420 s cintill ator and S10362-33-050 MPPCs. 

Figure [T2 . 2 1 shows the hit rate per tile as a func- 
tion of the position along the beam direction for 
phase II, including also tracks which pass several 
tiles, as well as background events e.g. from the 
collimators. These rates can be handled by the 
proposed readout electronics (see section 13 1. In 
order to further reduce the hit rate, and therefore 
the signal pileup, optimizations like adjusting the 
tile angle to match the mean incident angle of the 
tracks are currently studied. 

The time-stamps for the individual hits are as- 
signed using the signal waveforms generated by 
the SiPM simulation. A simple fixed threshold 
method is used which approximates the behaviour 
of the STiC chip. The same method can also be 
applied to the output of the DRS chip. 

The optimal threshold was determined to be 
~ 10% of the mean signal amplitude (« 350 
pixels). Due to the varying signal amplitudes and 
the fixed threshold, the individual time-stamps 
have to be corrected for the timewalk effect. Fig- 
ure 112.31 shows the timewalk as a function of the 
Time-over- Threshold . 




ToT [ns] 

Figure 12.3: Timewalk correction. 



Figure |12.4| shows the distribution of the time- 
stamps relative to the true hit time given by the 
simulation for phase II. A fit of a Gaussian func- 
tion yields a time resolution of a t = 45 ps. The 
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Figure 12.4: Timing for phase II. The timestamps 
have been corrected for timewalk. 



non-Gaussian tails of the distribution come from 
signal pileup. 

The efficiency for assigning a time-stamp to an 
isolated track within ±3ert of the true hit time 
is i=a 100%. Inefficiencies at the edge of a recurl 
station, where a track deposits too little energy in 
a tile to pass the threshold can be neglected. 

However, the efficiency is also influenced by sig- 
nal pileup. For the simple fixed threshold model 
used in this analysis, an overall efficiency of sa 
99.5% for phase I and « 98.0% for phase II is 
achieved. With a more sophisticated peak-finding 
algorithm, it should be possible to resolve pileup 
signals which occur within a time interval of a 
few nanoseconds and consequently achieve an ef- 
ficiency of w 100%. Also, using a SiPM with a 
fast output, like the SensL MicroFB series, will 
significantly reduce the pileup and thus increase 
the efficiency. 



12.3 Time Resolution 
ments 



Measure- 



First measurements of time resolution of a single 
lxlxl cm 3 tile and a Hamamatsu S10362-33-050 
MPPC have been done using an oscilloscope with 
a bandwidth of 1 GHz. The scintillator material 
used in this measurements is NE110. Compared 
to BC420, this scintillator has a lower light yield 
and a slightly slower rise and decay time; using 
a more suited scintillator will thus improve the 
results. The scintillation light was triggered using 
a pulsed UV laseiQ 



Figure |12.5| shows the measured time resolu- 
tion as a function of the signal amplitude for a 
threshold of 10%. For a typical signal amplitude 
of s» 350 pixels, a time resolution of sa 45 ps 
is achieved, which is consistent with the simula- 
tion results. Figures [12. 6| and |12.7| show the de- 
pendence on the applied threshold and SiPM bias 
voltage. 

These results show, that a time resolution well 
below 100 ps can be achieved. It is expected, the 
resolution will not be degraded using the DRS or 
STiC readout chip, which will be verified in future 
measurements . 
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Figure 12.5: Time resolution as a function of the 
signal amplitude for a threshold of 10%. The amp- 
litude expected in the tile detector is sa 500 pixels. 
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Figure 12.6: Time resolution as a function of the 
discrimination threshold for 500 pixel signal amp- 
litude. 

1 The scintillator response to the UV light is a good ap- 
proximation of the response to electrons. 
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12.4 Detector Prototype 

The design of a prototype of a tile detector mod- 



ule has been started (see Figure 12.1 1 and the 
construction is expected to be completed within 
2013. In this prototype the mechanical design of 
the support structures and readout infrastructure 
will be tested and measurements of the detector 
performance will be carried out. 



Figure 12.7: Time resolution as a function of the 
SiPM bias voltage for 500 pixel signal amplitude 
and a threshold of 10%. 
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Chapter 13 




Data Acquisition 



1116 Pixel Sensors 



up to 108 

800 Mbit/s links 



42 FPGAs 



gets to see the complete detector information for 
a select time slice. See Figure 13.2 for an over- 
view of the readout scheme and Figure |13.1| for 
the scheme at detector start-up. 



13.2 Occupancy 
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Figure 13.1: Mu3e readout scheme for the start-up 
detector. 



13.1 Overview 

The Mu3e data acquisition system works without 
a hardware trigger on a push basis, i.e. the de- 
tector elements continuously send hit information 
to the data acquisition (DAQ) system. The DAQ 
consists of three layers, namely front-end FPGAs, 
read-out boards and the filter farm. The topology 
of interconnects is built such that every farm PC 



The bandwidth requirements of the data acquis- 
ition are largely determined by the expected de- 
tector occupancy, as all the Mu3e sub-detectors 
produce zero-suppressed output. 

The occupancies shown are obtained with the 
full simulation running at a muon stop rate of 
2 • 10 9 Hz (2 • 10 8 Hz for phase I) and pessimistic- 
ally estimating the beam related background by 
loosing another 4 • 10 9 Hz (4 • 10 8 Hz) of muons 
along the beam line. Figures [13. 3| and |13.4| show 
the expected number of hits per 50 ns frame in 
the pixel detector. Figures [l 3 . 5 1 and 1 1 3 . 6| show the 
same for the fibre detector. The distribution of 
the occupancy over the pixel sensors is shown in 
Figures pTf] and [I3~8l 

13.3 Front-end 

13.3.1 Pixel detector 

The pixel sensors contain electronics for hit de- 
tection and time as well as address encoding. All 
hits assigned to the same (20 MHz) time-stamp 
constitute a frame. The sensors collect the data 
of 16 frames into a superframe and send it off chip 
via an 800 Mbit/s low- voltage differential signal- 
ing (LVDS) link. The signals travel over a max- 
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Figure 13.2: Overall Mu3e readout scheme 
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Figure 13.3: Number of pixel hits in the central 
detector per 50 ns frame in phase I running. 
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Figure 13.4: Number of pixel hits in the complete 
detector per 50 ns frame in phase II running. 
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Figure 13.5: Number of fibre hits per 50 ns frame 
in phase I running. 
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Figure 13.6: Number of fibre hits per 50 ns frame 
in phase II running. 




Layer 2 



2 4 



Layer 3 



2 4 6 8 10 12 14 16 



Layer 4 




Figure 13.7: Occupancy in 50 ns frames of the 
central pixel sensors for phase II running. The 
axes enumerate sensor numbers. For phase I, the 
occupancy numbers have to be scaled down by a 
factor 10-20. 



imum of 18 cm on a Kapton flex-print to the edge 
of the sensitive area, where they are amplified by 
a driver chip. The Kapton prints then connects to 
a PCB located between the recurl layers and the 
beam-pipe. On this PCB, up to 72 LVDS links 
are fed into a FPGA. The FPGA provides buffer- 
ing and collects a long stream of frames (at least 
1024) into a frametrain. The assembled data are 
then output to 8 3 Gbit/s links, such that the data 
of one frame-train are sent on two links. On the 
PCB, the signals are converted to optical and sent 
off-detector via fibres. An additional pair of op- 
tical links per FPGA is required for slow control 
and monitoring. 



Hardware 

The requirements for the on-detector FPGAs can 
be met by mid- or even low-price devices (such as 
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Tabic 13.1: Pixel readout requirements (Phase IB without recurl stations). 
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Table 13.2: Pixel readout requirements (Phase II), for the recurl stations only the sum of average hits 
per layer and FPGA— >RO capacity is given. 



the ALTERA Cyclone IV family or the XILINX 
Artix VII family) . The FPGAs are to be mounted 
on PCBs that are placed between the recurl layers 
and the beam-pipe. 

Firmware 

The main task of the on-detector FPGAs is col- 
lecting the relatively short time slices of 16 clock 
cycles assembled on the pixel chips to the long in- 
tervals treated by the individual filter farm PCs. 
During this buffering, the hits can be time ordered 
inside a slice and the protocol overhead can be re- 
duced. In addition, hits can be clustered. 

A further task for the first line of FPGAs is the 
configuration and monitoring of the pixel chips. 
A 32 bit histogram of the hit counts in a single 
sensor however requires 256 kB of memory, thus 
exceeding the capacity of the devices; an ex- 



ternal memory interface would significantly in- 
crease the pin count and the PCB complexity; the 
histograming task is thus deferred to the readout 
boards. 

These tasks arc all fairly standard and FPGAs 
that fulfill the bandwidth requirements for the in- 
and output channels do provide enough logic for 
implementing them. 



13.3.2 Timing detector 

For the timing detectors, three readout schemes 
are currently under investigation: One based on 
a further development of the DRS switched ca- 
pacitor array developed at PSI, one based on the 
STiC chip developed at KIP, Heidelberg Univer- 
sity and one with FPGA-based TDCs (mainly for 
the fibres). 
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Figure 13.8: Occupancy in 50 ns frames of the re- 
curl pixel sensors for phase II running. The axes 
enumerate sensor numbers. 



DRS sampling readout 

The readout of the tile and fibre detectors requires 
high rate capability and extremely good timing 
resolution. To achieve an overall detector tim- 
ing accuracy below 100 ps, the associated elec- 
tronics needs to be at least a factor of two bet- 
ter, i.e. 50 ps. The high rate environment causes 
significant pile-up, which limits the usage of con- 
ventional techniques such as discriminators and 
TDCs. Therefore one option to read out the tiles 
is with the well-established waveform digitizing 
technology developed at PSI, which is in use since 
many years in the MEG experiment. It is based 
on the DRS4 switched capacitor array, which is 
capable of sampling the SiPM signals with up to 5 
Giga samples per second (GSPS) with a resolution 
close to 12 bits. It has been shown in the MEG 
experiment that this technology allows a timing 
accuracy in the order of 40 ps across many thou- 
sand channels. The knowledge of the exact wave- 



form of an event is very well suited to detect and 
suppress pile-up. 

The tile readout electronics could be placed out- 
side the detector in special crates connected with 
a few meters of cable. This simplifies the design 
and maintainability, while not compromising the 
signal quality dramatically. 

A principal limitation arises from the DRS4 
chip, which is capable of only a limited event rate 
of about 100 kHz. While this will be sufficient 
for Phase IB, it has to be improved for Phase 
II. Therefore a new development has been star- 
ted to design a new version of this chip. The 
DRS5 chip will use an internal analog memory 
(FIFO) to work in a dead-time less fashion up to 
an event rate of about 5 MHz. A critical part of 
the DRS5 chip which is the inverter chain oper- 
ating the sampling circuitry has already been de- 
signed in the new 110 nm CMOS technology and 
submitted. First test results are expected begin- 
ning of 2013. The dead-time less operation of this 
chip will be combined with higher sampling speed 
(10 GSPS) and a better timing accuracy, allowing 
for a time measurement well below 10 ps. 

In order to limit the amount of data to be read 
out, the FPGA connected to the DRS chip will 
already analyze the waveform and extract its ma- 
jor parameters like time and amplitude. Only a 
prescaled subset of events will contain the full de- 
tector waveform in order to cross-check the ana- 
lysis algorithms in the FPGA. Methods have re- 
cently been published which obtain the timing in- 
formation by cross correlation or cubic interpol- 
ation with an accuracy of about 1/1 0th of the 
sampling interval, which would be 10 ps in the case 
of 5 GSPS. 



STiC readout 

The STiC chip offers an alternative to the DRS5 
readout. STiC is a mixed mode 16-channel ASIC 
chip in UMC 0.18 urn CMOS technology designed 
for SiPM readout with high time resolution. It 
is developed for Time-of-Flight measurements in 
high energy physics and medical imaging, in par- 
ticular the EndoTOFPET-US project. The chip 
has a differential structure, however, it supports 
both differential and single-ended connection of 
SiPMs. A 6 bit DAC allows to tunc the voltage 
at each input terminal within « 1 V. In this way 
the SiPM operating voltage can be adjusted and 
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makes programming much more fickle. It has how- 
ever been shown that resolutions of O(10 ps) can 

Whilst 



be achieved with this technology [149] 
definitely not the optimal solution, FPGA based 
timing could serve as a low-cost, low-risk solution 
for phase IB running. 




13.4 Read-out links 
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Figure 13.9: Dual threshold discrimination for en- 
ergy and timing information. 



temperature and device-to-device fluctuations can 
be compensated. 

The time and charge information of the signal 
are encrypted into two time-stamps which are ob- 
tained by discriminating the signal with two differ- 
ent thresholds (see Figure 13.91. The thresholds 
can be tuned in a range of w 0.2 — 15 pixel sig- 
nals for the timing and up to 200 pixel signals for 
the charge. The time-stamps are then processed 
by an embedded TDC module with a resolution 
of < 20 ps. A special linearization method is im- 
plemented to obtain a linear charge response in a 
very wide range. With the chip, a time-resolution 
of « 50 ps was measured for a 10 pixel signal of 
a MPPC S10362-33-50 without scintillator. For 
typical signal amplitudes of 0(100) pixels, which 
are expected in the tile detector, the timing jitter 
of the chip is negligible. 

The data rate of the current chip is limited to 
Rj 100 kHz per channel. However, the chip will be 
modified to allow for data rates of w 500 kHz per 
channel within 2013/2014. Until 2016, the data 
rate will be further increased to several MHz in 
order to match the requirements for phase II. 



FPGA based readout 

A further alternative for the fibre readout is the 
use of time to digital converters implemented in 
FPGAs. Resolutions of 0(1 ns) can be achieved 
fairly cheaply; much better performance requires 
the use of carry chain techniques, which greatly 
reduces the number of channels per FPGA and 



In total there are three different types of read-out 
links in the Mu3e data acquisition system. The 
same type of links can be used for a small num- 
ber of slow and fast control links in the opposite 
direction. 

The data from the MUPIX chips will be trans- 
mitted to the front-end FPGAs via LVDS links 
at 800Mbit/s, which is a quasi industry stand- 
ard. The fast serializers and LVDS drivers for the 
MUPIX chip will be adopted from a similar chip 
design by a group from Bonn 150 . The link will 



be physically implemented as a matched differen- 
tial pair of aluminum traces on the sensor flex- 
print, shielded by ground lines on both sides. It 
is foreseen to have a first set of LVDS repeaters 
just outside the acceptance of the detector. Very 
compact commercially available LVDS repeaters 
have already proven to be radiation tolerant. Out- 
side the acceptance, the flex-print will then be of 
multilayer type, with differential signal inner lay- 
ers, shielded by ground planes. The connection to 
the PCB housing the front-end FPGA is made by 
small form factor connectors, which support high 
bandwidth. Lastly the routing on the front-end 
FPGA PCB will be done with matched differen- 
tial copper traces going to the pre-defined LVDS 
inputs of the FPGA itself. 

There will be two types of optical high speed 
data links. The first one is going from the front- 
end FPGAs to the read-out boards, the second 
from the read-out boards to the FPGA PCIc 
boards in the event filter farm PCs. The optical 
links from the front-end FPGAs to the read out 
boards have a bandwidth of 3.125 Gbit/s, which 
fits well the FPGA specifications. Each FPGA has 
eight (or more) fast transceiver blocks. The op- 
tical interface can be implemented using radiation 
tolerant components developed by the Versatile 
Link group |151| . Especially the 12- way laser ar- 
ray with MTP/MPO standard connector would 
give a very compact form factor, the fall back solu- 
tion would be a radiation hard dual laser transmit- 
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ter system, the VTTx module developed at CERN 
152 . The data laser will have a wavelength of 



850 nm and the optical fibre is of 50/125 multi- 
mode type, since this is a standard both in in- 
dustry and in particle physics detector readout. 
In the case of the 12- way emitter the fibres will go 
to a splitter and patch panel in order to combine 
fibres going to the same read-out board. In the 
case of dual transmitters the fibres can be com- 
bined in groups of 12 using break-out fibre cables. 
The fibres traveling from the experimental area 
to the filter farm counting house can be further 
combined in cables containing 8 x 12 fibres going 
to the same sub-farm. An optical patch panel for 
each sub-farm will allow to connect the long dis- 
tance 96-fibre cables to 12 fibre patch cords going 
to the readout boards. The readout boards have 
2 or 3 12-way optical inputs, which are populated 
with commercial receiver modules. The high-end 
FPGAs used on the read-out card have up to 66 
high speed transceiver blocks, which will be used 
for the de-serialization of the data. 

The last type of link connects the read-out cards 
to the FPGA PCIe boards in the event filter farm 
PCs. This optical link will be implemented as 
lOGbit/s high speed link. Since each read-out 
board is connected to every sub-farm PC with one 
of the high speed links, they are point to point 
single fibre links. The fibres are of 50/125 multi- 
mode type operated at 850 nm to stay compatible 
with the links from the detector. If 12-way optical 
transmitters and receivers running at lOGbit/s 
can be purchased, they will be used for this high 
speed link. Otherwise single fibre transmitters 
and receivers will be used both on the read-out 
board and the FPGA PCIe board side. In prac- 
tice the high speed optical receivers for the FPGA 
PCIe boards could be on optical mezzanines such 
as the 8 channel St. Luce card developed by TU 
Dortmund 



153 



13.5 Read-out cards 

The main tasks of the read-out cards is to act as 
switches between the front end and the on-line re- 
construction farm and to act as buffers between 
the synchronous front end and the asynchronous 
back end. The board design and choice of FP- 
GAs is dominated by the number of fast links re- 
quired. We plan to adapt an existing develop- 



compute nodes 155 , which would both fulfill our 
needs. 



13.6 Event filter interface 

The filter farm PCs will be equipped with FPGA 
cards in PCIe slots and optical receiver daughter 



cards, as described in more detail in section 14.2 



13.7 Data collection 

The filter farm will output selected events at a 
data rate in the order of 50 MBytes/s in total. 
This data rate is low enough to be collected by a 
single PC connected to the filter farm by com- 
mon GBit Ethernet and written to local disks. 
Then the data will be transferred to the central 
PSI computing center, where it is stored and ana- 
lyzed. For the central data acquisition the well 
established MIDAS (Maximum Integrated Data 
Acquisition System) [156] software package will 
be used. This software is currently used in several 
major experiments such as the T2K ND280 de- 
tector in Japan |157| , A LPH A at CERN and the 
MEG experiment at PSI 158 . It can easily handle 



ment, e.g. LHCb TELL1 cards 154 or PANDA 



the required data rate, and contains all necessary 
tools such as event building, a slow control system 
including a history database and an alarm system. 
A web interface allows controlling and monitoring 
the experiment through the Internet. The farm 
PCs will use MIDAS library calls to ship the data 
to the central DAQ PC. The framework also offers 
facilities to send configuration parameters from a 
central database (the "Online DataBase" or ODB) 
to all connected farm PCs and to coordinate com- 
mon starts and stops of acquisition (run control) . 

For the purpose of monitoring and data qual- 
ity control of the experiment the MIDAS system 
offers taps to the data stream for connections of 
analysis and graphical display programs. The out- 
put of graphical user interface programs can be fed 
back into the web interface of the MIDAS system 
so the experiment can be monitored also remotely 
with just a Web browser. 

13.8 Slow control 

The slow control system deals with all "slow" data 
such as high voltages for the SiPMs and silicon 
sensors, ambient temperatures and pressures. For 
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Figure 13.10: SCS-2001 unit as part of the 
MSCB slow control system. This unit has 64 in- 
put/output channels, which can be configured via 
plug-in boards as digital or analog channels. Many 
plug-in boards exist already such as PT100 tem- 
perature sensor readout cards, analog high resol- 
ution inputs (24 bit resolution) , valve control out- 
puts and many more. 

the configuration and control of the silicon pixel 
sensors the JTAG standard [132] will be used. It 
is planned to use the MIDAS Slow Control Bus 
(MSCB) system |159| to link all distributed con- 
trol and monitoring devices into a single system. 
The MSCB system is also well established at sev- 
eral laboratories. It uses a serial differential bus 
for communication, and simple micro controllers 
in all control devices. The micro controllers per- 
form local control loops such as high voltage sta- 
bilization, and send measured values to the cent- 
ral DAQ system for monitoring. Many devices 
already exist for this system, such as the SCS- 
2001 unit shown in |13.10| Since the system was 



developed at PSI, it can be quickly adapted to 
new hardware. The high voltage control for the 
SiPMs can for example be directly integrated into 
the carrier boards holding the SiPMs, thus elim- 
inating the need for high voltage cables. The op- 
timized protocol of the MSCB system allows the 
monitoring of many thousand channels with re- 
petition rates in the 100 ms range, which will be 
more than enough for this experiment. 

In addition to the MSCB system, the MIDAS 
slow control package contains interfaces to the PSI 
beamline elements via the EPICS system 160 



This allows monitoring and control of the beam- 
line from the main DAQ system, which has been 
proven very versatile in other experiments using 
this scheme. 

All slow control data will be stored in the his- 
tory system of the MIDAS system, so that long 
term stabilities of the experiment can be effect- 
ively verified. The slow control data is also fed 
into the main event data stream, so that any off- 
line analysis of the event data has this data avail- 
able. 

A special case is the configuration of the pixel 
detectors, which require many million parameters, 
like the trim-DAC values for each pixel. Since the 
amount of data here is considerably larger than for 
all other systems, an extension of the slow control 
system is planned. A dedicated program manages, 
visualizes and exchanges the pixel detector con- 
figuration parameters between an optimized data- 
base and the pixel hardware. In this way the time 
required to configure the pixel detectors can be 
minimized, while this program is still connected 
to the main DAQ system. It can be synchronized 
with run starts and stops, and can inject pixel 
monitoring data periodically into the event data 
stream for offline analysis. 
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Chapter 14 




Online Event Selection 



14.1 Selection Algorithms 

As in the final analysis, event selection in the filter 
farm can rely on the coincidence of three tracks in 
time and vertex and on their kinematics. Espe- 
cially for high rate running, coincidence in time 
in the fibre detector is not sufficient to reduce the 
data rate by three to four orders of magnitude. 
Thus a track reconstruction will be required. The 
triplet based multiple scattering fit described in 
chapter [16] is well suited for online implementa- 
tion and current GPUs can perform 10 9 triplet 
fits per seconcQ thus already fulfilling the needs 
of Mu3e up to at least medium intensity (few 10 8 
muons/s) running. 

Triplets of the tracks thus reconstructed can 
then be fit to a common vertex. Even loose ver- 
tex requirements can give a 10 3 reduction factor 
at 2-10 9 Hz muon rate and 10 4 - 1 5 f or th e 
phase I experiment (see Figures 14.1 and 14.21. 



Combining the vertexing with modest kinematic 
requirements (e.g. on the three-particle invariant 
mass or the planarity) should produce the required 
data reduction, leaving the timing information as 
a valuable offline cross-check (and obviating the 
need for online timing reconstruction). 



14.2 Hardware Implementation 

The data will arrive on the farm PCs via op- 
tical links on a PCIe FPGA board. The FPGA 
will perform the event building and buffering and 



also allows to run simple clustering and sorting 
algorithms. The event data are then transferred 
via DMA over the PCIe 3 bus^to the memory of 
a graphics processing unit (GPU), where the se- 
lection algorithms are run. The GPU then posts 
selected events and monitoring data to the main 
memory of the PC, from where the CPU ships 
it via Ethernet to the central data acquisition 
computer running the MIDAS software. At that 
computer, the data streams from the farm PCs 
are combined into a single data stream, combined 
with various slow control data, compressed and 
stored. 

For the receiver FPGA cards, evaluation boards 
16il, or ALTERA (Fig- 



from either XILINX 
ure 



14.3 1 162|163 or similar hardware built by the 



1 As tested on a AMD Radeon 6990 using OpenCL under 
Linux. 



collaboration could be used in conjunction with 
daughter boards with the optical receivers (sim- 
ilar to e.g. the optical receiver boards used in the 
LHCb readout electronics |164| ). The maximum 
data rate over the PCIe 3.0 bus is 16Gbyte/s, 
amply sufficient for phase ^] For the full phase II 
rate, the raw link speed is still sufficient, would 
however have to be fully and efficiently used. The 
PCIe 4.0 standard, doubling this rate, should be- 
come commercially available around 2017, com- 
patible with phase II running; alternatively, the 
number of farm PCs could be increased. 



2 Note that PCIe is actually not a bus protocol, but 
offers switched point-to-point connections. The bus desig- 
nation is due to the software-side backwards compatibility 
to the original PCI bus interface. 

3 For phase I running, the FPGA-GPU link can also be 
implemented on PCIe 2.0 (max. 8 Gbyte/s), which is better 
supported on currently available FPGAs. 
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Figure 14.2: Fraction of 50 ns frames containing 
three vertices consistent with e + e~e + inside a 
given distance for a muon stop rate of 2 • 10 9 Hz for 
680'000 simulated frames. In the top plot, every 
crossing of a simulated electron/positron track is 
counted as a vertex; charge assignments are made 
purely on the apparent curvature, i.e. recurling 
positrons are counted as electrons. In the bottom 
plot, only true electrons are counted. 




Figure 14.3: ALTERA Stratix IV PCIe develop- 
ment board. 
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Figure 14.1: Fraction of 50 ns frames containing 
three vertices consistent with e + e~e + inside a 
given distance for a muon stop rate of 2 ■ 10 8 Hz 
for 3.15 million simulated frames. In the top plot, 
every crossing of a simulated electron/positron 
track is counted as a vertex; charge assignments 
are made purely on the apparent curvature, i.e. re- 
curling positrons are counted as electrons. In the 
bottom plot, only true electrons are counted. 

The GPU boards will be obtained commercially 
as late as possible in order to profit from the fast 
developments and sinking prices. As far as raw 
floating point throughput is concerned, current 
high-end GPUs already pack enough power for 
high rate running 165[|166 . Newer cards are how- 
ever expected to offer higher memory bandwidth 
and better caching. Also the performance of the 
driver software (especially as far as the PCIe 3 bus 
is concerned) and the GPU compilers is expec- 
ted to improve. The two GPU vendors AMD and 
NVIDIA offer fairly different architectures; which 
one performs better depends a lot on the details of 
the algorithm to be implemented; we arc currently 
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performing tests with both architectures and will 
choose a vendor once we have a mature imple- 
mentation. 

We currently plan to host the farm PCs in in- 
dividual tower casings, ensuring enough space for 



the FPGA board and the high end GPU whilst 
allowing for air cooling. At load, each tower will 
consume around 0.5 KW, so adequate cooling of 
the counting house is essential. 
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Chapter 15 




Simulation 



This chapter describes the Geant4 124 167 



based simulation used to derive the figures and 
plots in this proposal. 



15.1 Detector geometry 

15.1.1 Beam delivery 

In the simulation, the beam is started 3 m in 
front of the target inside a beam transport solen- 
oid. Beam particles are generated with a profile 
and momentum spectrum like the one observed in 
MEG. 1.5 m before the target, the beam enters the 
main solenoid and shortly before the target it exits 
the beam vacuum through a thin window. Along 
the beamline, two thick lead collimators reduce 
the beam to the target size. For an overview of 
the simulated beamline elements, sec Figure [l~5T2| 
In this simple setup, about a third of the gen- 
erated muons decay in the target, which, whilst 




Figure 15.1: 
detector. 



Wire frame view of the simulated 



not very efficient, gives a conservative estimate of 
beam-induced backgrounds. 

15.1.2 Target 

The target is simulated as a hollow aluminium 
double cone supported by three nylon strings at 
each end and a nylon string along its axis, see also 
chapter [9] 

15.1.3 Pixel detector 

The pixel detector is simulated as 50 um of sil- 
icon on top of 15 um of aluminium representing the 
traces on the flexprint (covering half the available 
area) on top of 50 um of Kapton, with the silicon 
offset such that an overlap with next sensor is cre- 
ated, see Figure [T5~i"| Half a millimeter of the pixel 
sensor at the edge is assumed to be inactive, the 
rest is divided into 80 x 80 um 2 pixels. The sim- 
ulated sensor layers are supported at their ends 
by plastic and aluminium structures modeled on 
those in the mechanical prototype shown in Fig- 
ure nTT2lH 

15.1.4 Scintillating fibres 

The fibre detector is simulated as three circular 
layers of 250 um scintillating fibres in the main 
simulation. A detailed simulation including op- 
tical transport and the effect of fibre cladding and 
coating also exists, see section [1 1. 3 1 The results of 
the detailed simulation regarding light yield and 
propagation times will eventually be fed back into 
the main simulation in a parameterized form. The 
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Mu3e simulation beamline 




Transport Solenoid B, = 1 .5T Field only (no material in Geant) 




z = -3000 mm 
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Figure 15.2: Beamline in the current simulation. 



100 mm 



Figure 15.3: Geometry of the detector in the simulation. The top half only shows active (sensor) 
volumina, whereas the bottom half only shows support structures. 
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Figure 15.4: Pixel detector simulation geometry 
for the innermost layer. The sensor is shown 
in black, the aluminium traces in blue and the 
Kapton support in orange. Note that all thick- 
nesses are stretched by a factor of 2. 



Geant4; in a test setup they were compared to the 
single scattering model, see Figure [T5. 5 1 The best 
overall description is obtained from the Urban- 
Model 168 at large step widths, which also has 



the shortest computation times. In the helium gas 
on the other hand, none of the parameterizations 



performs adequately, see Figure 15.6 



We plan to verify the simulation results with 
beam telescope measurements in 2013, which 
should also lead to a usable parameterization of 
multiple scattering in gases. 



response of the silicon photomultipliers is simu- 
lated by the GosSiP package [1481. In the simu- 



lation, the fibres are supported at both ends by 
massive aluminium rings. 

15.1.5 Tile detector 

The simulated tile detector consists of plastic scin- 
tillator tiles mounted on an aluminium tube. Also 
here, a separate detailed simulation including light 
transport and silicon photomultiplier response is 
available and will have to be fed back into the 
main simulation in a parameterized form. 



15.3.2 MUON DECAYS 



Michel decay 



Geant4 implements the Michel decay including po- 
larization of both the muon and the positron based 
and 



on 



169 



170 



The spectra of the neutrinos 
do not follow the physical distribution, this does 
however not affect the simulation for Mu3e. Some- 
what more worrying is the fact that the Michel 
matrix element contains radiative corrections but 
is not clearly separated from the radiative decay 
matrix clement. 



15.2 Magnetic field 



Radiative decay 



The simulated magnetic field can be read from 
arbitrary field maps or generated in the code via 
integration over current loops. The propagation 
of muons in the field includes spin tracking. For 
the simulations shown in this report, the field was 
generated from 100 current loops spaced equally 
over 3 m, with currents normalized such that the 
longitudinal component of the field in the center 
of the target is 1 T, supplemented by a 1.5 T field 
in the center of the beam transport solenoid, see 



section 15.1.1 and Figure 15.2 



15.3 Physics Processes 

15.3.1 Multiple Coulomb scattering 

Multiple coulomb scattering is the main limiting 
factor for the resolution of the experiment; an ac- 
curate simulation is thus crucial. The best res- 
ults are obtained by simulating each individual 
scattering, which however results in prohibitively 
large computing times. A large selection of mul- 
tiple scattering parameterizations are available in 



The radiative decay of the muon was implemented 
in Geant4 by the TWIST collaboration [m] based 



on 



172 



The code does not describe the neutrino 
spectra and avoids the collinear and infrared sin- 
gularities by sampling the matrix element assum- 
ing a finite integral. 



Radiative decay with internal conversion 

The radiative decay with internal conversion is 
simulated using the hit and miss technique on 
events generated evenly in phase space using the 



RAMBO c ode [173 and applying the matrix ele- 
ment from 116 . Unfortunately, there is currently 
no polarized version of the matrix element avail- 
able and thus the simulation is unpolarized. The 
hit and miss technique is very expensive in terms 
of computation time, if the complete phase space 
is to be simulated (as the matrix elements varies 
by more than 16 orders of magnitude), this can 
however be overcome by restricting the simulation 
to regions of particular interest, e.g. high invariant 
masses of the visible particles. 
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Figure 15.5: Comparison of multiple coulomb scattering models in different scattering angle ranges. The 
scatterer is a single silicon-Kapton assembly shot at at a right angle with 30MeV positrons. The black 
dots and the green line show the single scattering model which serves as a reference; as expected, the 
single scattering model is not affected by the Geant step size. Of all the parameterizations, the Urban 
model with a step size that treats each bit of material as a single volume performs best. 





0.15 0.2 0.25 0.3 0.35 0.4 0.45 




Figure 15.6: Comparison of multiple coulomb scattering models in m of helium gas for different scattering 
angle ranges. The test particles are 30 MeV positrons. The black dots and the green line show the single 
scattering model which serves as a reference. All the parameterizations are unfortunately inadequate. 
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Figure 15.7: View of a simulated frame. 



Signal 

The signal kinematics are highly model- 
dependent, see chapter [4] If not otherwise 
noted, we have used three particle phase space 
distributions in the simulation, following the 
practice of SINDRUM and earlier experiments. 

Special decays 

The simulation allows the simulation of overlap 
decays, where we force more than one muon decay 
to happen at a single vertex. Thus we can simu- 
late the accidental backgrounds arising e.g. from 
the overlap of an internal conversion decay and a 
Michel decay without having to generate in excess 
of 10 16 frames. 



15.4 Time structure 

As the Mu3e experiment operates with a quasi 
continuous beam, the paradigms of bunch cross- 
ing and event from collider experiments do not 
apply; they have however informed the design of 
the Geant4 package. In our simulation, particles 
are assigned a 64 bit ID, which is unique over 
a run and thus conserves identities and mother- 
daughter relationships across the boundaries of 
read-out time frames. Before each step of track- 
ing the particle through the detector, it is checked 
whether the particle has just crossed into the next 
time slice. If so, its information is stored, its 
time relative to the time slice adjusted and track- 
ing deferred to the next slice. Thus we ensure 
that we correctly treat muons stuck in the tar- 
get before decaying and decay products crossing 
read-out frame boundaries while traversing the de- 
tector. In order to simulate a steady state, where 
approximately the same number of muons enter 
the target and decay, the first 5 ms of simulation 
running, during which the target is loaded, are 
usually thrown away and not used in occupancy 
or efficiency studies. 

Currently not simulated are effects of the 
40 MHz structure of the primary proton beam on 
the time structure seen in the detector; if this 
would be needed, a measured structure could eas- 
ily be superimposed on the generation of muons 
in the simulation framework. 
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Chapter 16 




Reconstruction 



16.1 Track Reconstruction in the 
Pixel Tracker 

A precise track reconstruction of electrons is of 
highest importance for the identification of the 
jU — > eee decay with a sensitivity of 1 out of 
10 15 (10 16 ) ordinary Michel decays in phase I (II), 
which have to be suppressed by 16 orders of mag- 
nitude. 

Due to the high rate and the resulting high oc- 
cupancy especially at phase II of the project with 
up to 100 tracks per readout frame, the recon- 
struction algorithm has to deal effectively with 
the combinatorial background in order to reduce 
the fake rate, i.e. the rate of wrongly reconstruc- 
ted tracks, to an acceptable level. The combin- 
atorial problem is not only due to the high rate 
but also due to the large bending of the low mo- 
mentum electrons in the strong magnetic field of 
B = 1 T, which, depending on the position and 
flight direction can make several turns in the de- 
tector (recurlers) . Hit combinations can span over 
distances of more than half a meter. Hits of re- 
curling tracks are found on opposite sides of the 
detector and still have to be correctly combined by 
the reconstruction program. This is of particular 
importance for the determination of the flight dir- 
ection and therefore charge of the particle. Only 
for a fully reconstructed track the time informa- 
tion provided by the time of flight system can be 
correctly applied. 

As the full detector readout is triggerless, 
all muon decays have to be fully reconstructed 
already on filter farm level, setting high demands 



on the speed of the online track reconstruction al- 
gorithm. A further complication comes from the 
fact that the track resolution is dominated by mul- 
tiple scattering in the silicon pixel sensors and not 
by the pixel size, in contrast to most other exper- 
iments. Therefore, standard non-iterative circle 



fits of tracks 174 as used in high energy experi- 



ments can not be used here. 

In order to reduce multiple scattering, the 
number of sensor layers are reduced to a min- 
imum in the detector design which, unfortunately, 
also reduces redundancy for track reconstruction. 
Therefore, the track reconstruction also has to 
cope with a minimum of information provided by 
only four sensor layers. 



16.2 Track Fitting and Linking 

For the track reconstruction two different ap- 
proaches are followed in the Mu3e experiment, the 
broken line fit [175 176 and the 



177 



fast linear fit 
The broken 



based on multiple scattering 
line fit determines hit positions and scattering 
angles simultaneously and was implemented in 2D 
175) [178] and recently also in 3D [l76j[l79]. It 
is based on linearisation of a previous circle fit, 
works non-iteratively and provides the correlation 
between all fit quantities. The broken line fit, 
however, requires knowledge of the assignments of 
hits to tracks from a previous linking step. There- 
fore, the broken line fit can only be used in the 
final step of the track reconstruction, also because 
a previous track fit is required for the linearisation 
procedure. 
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The fast three-dimensional multiple scattering 




MS 



Figure 16.1: Sketch of the variables used in the 
multiple scattering fit. 
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Track polar angle 

Figure 16.2: Track types versus track polar angle 
for Michel decays in phase IA. 



(MS) fit 177 is based on fitting the multiple scat- 



tering angles at the middle hit position in a hit 
triplet combination, see Figure |16.1| In this fit, 
spatial uncertainties of the hit positions are ig- 
nored. This is a very good approximation for the 
Mu3e experiment as the pixel resolution uncer- 
tainty given by cr p i X ei = 80/ ^/l2 um is much smal- 
ler than the uncertainty from multiple scattering 
in the corresponding sensor layer. The MS-fit re- 
quires a detailed knowledge of the material distri- 
bution in the detector for the calculation of the 
scattering angle uncertainty. It minimises the azi- 
muthal and polar scattering angles at the sensor 
corresponding to the middle hit and exploits en- 
ergy conservatiorj^] The hit triplet trajectory, rep- 
resented by two connected helical curves, is de- 
scribed by the following two equations: 
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(16.2) 



The quantities, also shown in the sketch of Fig- 
ure |16.1[ are the following: i? 3 rj is the three di- 
mensional track radius, which can be directly re- 
lated to the momentum of the particle for a given 
magnetic field; <I>i ($2) are the bending angles of 
the first (second) arc and doi (^12) and zoi (^12) 
are the distances between the hits in the plane 
transverse and longitudinal to the solenoidal mag- 
netic field between the first and second (second 



x the energy loss in the Mu3e experiment is only about 
80 keV per sensor layer and can be neglected for track find- 
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Figure 16.3: Track types versus track polar angle 
for Michel decays in phase IB. 
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Figure 16.4: Track types versus track polar angle 
for Michel decays in phase II. 
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Figure 16.5: Tracking performance for Michel decays in phase IA. 



and third) hit, respectively. These equations can 
be linearised and solved in a fast non-iterative pro- 
cedure 11771. 



This linearized MS-fit is used as basis for the 
full reconstruction of tracks in the pixel detector. 
Tracks with more than three hits are fitted by sub- 
sequently combining several hit triplets. In the 

long tracks 



180 



current reconstruction program 
combining hits from several pixel layers are recon- 
structed first, then shorter tracks with fewer hit 
assignments are reconstructed. This procedure is 
repeated until no hits are left or no further tracks 
are found. Tracks with unresolved hit ambiguities 
are ignored in the following study to ensure high 
quality tracks with low fake rate. Also tracks with 
less than four hits combined are ignored. 

The number of hits linked to tracks depends on 
the single hit efficiency, which is assumed to be 
98% in the following studies, on the track direction 
(polar angle), the geometry of the pixel tracker 
and the geometrical acceptance of the detector, 



which is largest in phase II with four recurl sta- 
tions. The multiplicity of linked tracks as function 
of the polar angle is shown in Figures |16.2| |16.3| 
and |16.4| for the three different detector configura- 
tions. The multiplicity of linked hits is highest in 
the central region of the detector, where recurling 
tracks can be fully reconstructed. In the current 
version of the reconstruction program, only single 
turns of tracks are reconstructed, corresponding 
to a maximum of eight hits - four hits assigned to 
the outgoing trajectory and four hits assigned to 
the returning trajectory. In cases where one hit is 
missing, only seven hits are linked. For upstream 
or downstream going tracks the recurlers are not 
fully reconstructed, leading to tracks with mainly 
four linked hits or six linked hits depending on 
the detector setup and the number of implemen- 
ted recurl stations. It should be noted here that 
there is a big qualitative difference between the 
four-hit tracks and tracks with more than 4 hits 
(recurlers) . The bending radius of recurling tracks 
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Figure 16.6: Tracking performance for Michel decays in phase IB. 



can be measured with much higher precision due 
to the larger bending in the magnetic field. This 
leads to very different momentum resolutions for 
the different detector setups. 

In phase IA without Recurl Stations most tracks 
are reconstructed with only four linked hits. This 
yields a momentum resolution for reconstructed 
Michel electrons of about RMS(p) = 0.73MeV/c, 
see Figure |16.5| This resolution improves consid- 
erably to RMS(p) = 0.44MeV/c by addin g two 
inner Recurl Stations in phase IB (Figure 16.61 
and it improves further to RMS(p) = 0.28MeV/c 
by adding two inner and two outer Recurl Sta- 
tions in phase II Figure [16.7 1. The reconstruction 



of the recurling tracks, which provide high resol- 
ution momentum information, is crucial for the 
success of the experiment. 



16.3 Vertex Fitting 

A vertex fit will be used to precisely reconstruct 
the position of the muon decay for signal events. 
The vertex fit intrinsically checks the consistency 
of the assumption that all three signal candidate 
tracks originate from the same vertex. A common 
vertex fit of the three candidate electrons allows 
a suppression of the combinatorial background of 
Michel decays by a factor of about 3 • 10 s due to 
the very good pointing resolution of the pixel de- 
tector, which is mainly given by the multiple scat- 
tering at the first sensor layer. Instead of using 
a common vertex fit, also the distance between 
target impact points of reconstructed and extra- 
polated tracks can be used to discriminate signal 
against background. For simplicity the latter ap- 
proach is used here for the following sensitivity 
study. 
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Rec. Momentum - Gen. Momentum [MeV/c] Reconstructed track polar angle 

Figure 16.7: Tracking performance for Michel decays in phase II. 
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Chapter 17 




Sensitivity Study 




Figure 17.1: Hits in a simulated event at 2 ■ 10 9 Hz 
muon rate, viewed in the plane transverse to the 
beam. 

A sensitivity study based on a full simulation 
of the pixel detector is performed to estimate the 
background contribution at the different phases of 
the Mu3e experiment. The simulation includes a 
first material description of the detector, which 
varies between the different detector phases. The 
increasing level of detector instrumentation at 
later stages of the experiment will provide more 
information for particle tracking on one hand. 
On the other hand, the scintillating fibre tracker, 
which will be added at a later stage of the exper- 
iment when running at high muon rates, will also 



add multiple scattering and consequently increase 
confusion in the track linking step. In addition, 
the factor 10 and 100 higher increased particle 
rates at the Phase IB and II, respectively, will sig- 
nificantly increase the combinatorial background 
from Michel decays and other accidentals. A de- 
tailed simulation is therefore required to estimated 
the sensitivity reach of the different experimental 
phases. The simulated setups are summarized in 
Table [171] 

17.1 Simulation and Reconstruc- 
tion Software 

In order to estimate the maximum achievable 
sensitivities at the different stages of the project 
the track reconstruction program 1 180 is inter- 
faced to a) a full GEANT simulation of the de- 
tector and b) to a fast simulation program, which 
is mainly used for the following sensitivity studies. 
Both simulations were cross checked to give com- 
parable values for the track parameter resolution. 

In a first test it was proven that the track re- 
construction program is relatively insensitive with 
regards to changes of the single hit efficiency. Us- 
ing tight reconstruction criteria it is found that 
a 5 % loss of signal hits corresponds to an about 
5 % loss of tracks and an about 10 % loss of signal 
events. 

Noise hits also affect the track reconstruction. 
From noise studies of the MUPIX2 prototype chip 
it is known that the noise rate is very small (no 
noise hits were observed in an overnight run at 
a low threshold). To simulate the effect of noise 
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Phase IA 


Phase IB 


Phase II 


Recurl stations 





2 


4 


Scintillating fibres 


No 


Yes 


Yes 


Charge measurement 


No 


Yes 


Yes 


Target stopping rate 1 (Hz) 


2 • 10 7 


2 • 10 8 


2 • 10 9 



Table 17.1: Mu3e configurations considered for the sensitivity studies. Note that the detector acceptance 
is about. 50 % for signal events, given by the geometry of the central pixel detector. 
1 Highest possible rates are assumed in the simulation in order to be conservative. 



100 additional hits are randomly generated in the 
pixel detectors, corresponding to a noise rate of 
10~ 6 Hz per pixel. This rate will be remeasured 
as soon as a first large scale detector including a 
full readout chain becomes available. 

17.2 Signal Acceptance 

Signal events are generated by using a constant 
matrix element in phase space. This choice is 
motivated by LFV models with effective four- 
fermion contact interactions, which predict a con- 
stant matrix element. The so generated signal 
events are used for the following efficiency stud- 
ies. 

The resulting signal efficiency is mainly given 
by the acceptance of the experiment, which is in 
total about 50 % for the innermost detector lay- 
ers being placed at radii of approximately 2.0 cm 
and 3.0 cm. The total acceptance consists of two 
contributions, one related to the orientation of the 
decay plane with respect to the instrumented re- 
gions of the experiment and one from the min- 
imum bending radius (momentum), which can be 
reconstructed in the experimental setup. 

By moving the inner pixel layers closer to the 
target, the total acceptance can be slightly in- 
creased to about 55%. Further improvements of 
the total acceptance to about 60 % or more are 
only possible by using a longer inner barrel design 
and by reconstructing particles with a lower mo- 
mentum threshold. The latter can be achieved in 
two ways - a) by including highly bent trajectories 
in the track reconstruction or b) by decreasing the 
strength of the magnetic field, which however also 
compromises the momentum resolution. 

In the following, a conservative design is chosen 
with inner sensor layers placed at 2.0 and 3.0 cm 
and outer sensor layers placed at 7.0 and 8.0 cm, 
providing a signal acceptance of about 50 %. 



The resulting mass resolutions calculated from 
the reconstructed track parameters are shown in 
Figures 17.2||17.4 for the different phases. Com- 
pared to phase IA, the mass resolutions improve 
considerably by adding the recurl stations at 
phase IB and phase II and yield - despite the addi- 
tional material from the scintillating fibre tracker 
- resolutions well below lMeV. Of special import- 
ance for the background suppression is the strong 
reduction of the tails in the phase II setup. 



17.3 Selection 



17.3.1 Kinematic 
Events 



Selection of Signal 



The mass resolution can be further improved by 
rejecting badly reconstructed signal events, where 
one or several tracks are either wrongly recon- 
structed or suffer from large multiple scattering. 
A wrong measurement of the track momentum or 
direction affects momentum balance and leads to 
a measurable missing momentum. In the follow- 
ing we follow the strategy of the SINDRUM ex- 



periment 15 and define an acoplanar momentum 



vector Pacopi, which is obtained by projecting the 
vectorial sum of all track momenta into the decay 
plane defined by the three tracks. The correlation 
between the reconstructed invariant mass and the 
variable p a copi is shown in Figures 17.5 In ad- 



dition to the main spot at small values of p a copi, 
which originates from well measured signal events, 
two diagonal sidebands originating from wrongly 
measured signal events are visible. By applying 
the cut Pacopi < 1.4MeV/c, which will also be used 
to reject background events in the following, most 
of the wrongly measured signal events are rejec- 
ted. 

The resulting mass resolution plots are shown in 
Figure 17. 6f 17.81 which show by about 20 % im- 
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proved resolutions and significantly reduced tails. 
The expected mass resolutions and signal efficien- 
cies for the three stages of the experiment are sum- 
marized in Table 11731 

The main criterion for the kinematic selection of 
signal events is the reconstructed invariant mass. 
In a final analysis the number of signal and back- 
ground events will be determined by a fit of the 
invariant mass distribution, see also Figures [l7.12f 
|17.14| or by exploiting multivariate methods which 
include several estimators. For sake of simplicity 
the number of signal events is here determined 
from a 2-sigma mass window around the nominal 
muon mass. 



17.3.2 Reduction of the /x 
ground 



Back- 



The acoplanar momentum cut is also very effect- 
ive in reducing the dominating background from 
radiative events with internal conversion fi — > 
eeevv. These events are characterized by miss- 
ing energy carried away by the two undetected 
neutrinos. Most of these background events have 
either a small value of the reconstructed three elec- 
tron mass m eee or show some momentum imbal- 
ance. In particular the class of dangerous back- 
ground events, which have only little missing en- 
ergy, carried away by the two neutrinos, and are 
wrongly reconstructed such that the reconstruc- 
ted invariant mass matches the muon mass, show 
some significant momentum imbalance. Most of 
these background events are rejected by the cut 
Pacopi < 1.4MeV/c. 

The separation of signal and background in the 
plane p a copi versus m eee is shown for the three 
stages of the experiments in Figures 17.9p7.lT 



clear separation of signal and background events 
at a level of 1CP 17 — 1CP 16 muon decays is vis- 
ible. This separation improves with the upgraded 
experiment at phase IB and II. 

The projected invariant mass distribution of sig- 
nal and background after applying the acoplanar 
momentum cut is shown in the Figures |17.12f 
117.141 Theses distributions are used as basis for 
the estimated sensitivity calculation. 

17.3.3 Reduction of Accidental Back- 
ground 

The dominant contribution to the accidental back- 
ground comes from combinatorial background of 
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Figure 17.2: Reconstructed mass resolution for 
signal events in the phase IA configuration. 
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Figure 17.3: Reconstructed mass resolution for 
signal events in the phase IB configuration. 
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Figure 17.4: Reconstructed mass resolution for 
signal events in the phase II configuration. 
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Phase IA 


Phase IB 


Phase 11 


Michel decays: 
Efficiency (unpolarized) 
Momentum RMS 
Wrong charge fraction 


50.0% 
0.73MeV/c 
1.14% 


53.4% 
0.44MeV/c 
0.45% 


52.4% 
0.28MeV/c 
0.45% 


Signal: 
Reconstruction efficiency 
Energy sum RMS (reconstructed) 
Efficiency after selection 
Energy sum RMS (selected) 


39% 
1.12MeV/c 2 

26% 
0.91MeV/c 2 


46% 
0.65MeV/c 2 

39% 
0.47MeV/c 2 


48% 
0.52MeV/c 2 

38% 
0.42MeV/c 2 


Track d ca resolution (a) 


190 urn 


185 urn 


185 pm 



Table 17.2: Efficiencies and resolutions used in the sensitivity study. d ca is the distance of closest 
approach of a track to the beam line. The drop in the efficiency after selection for phase II is due to the 
larger combinatorial background. 
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Figure 17.5: Reconstruced mass versus acoplanar 
momentum for the phase II detector. 
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Figure 17.7: Reconstructed mass resolution for 
signal events after kinematic cuts in the phase IB 
configuration. 
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Figure 17.6: Reconstructed mass resolution for 
signal events after kinematic cuts in the phase IA 
configuration. 
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Figure 17.8: Reconstructed mass resolution for 
signal events after kinematic cuts in the phase II 
configuration. 
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Reconstructed Mass [MeV/c 2 ] 



Figure 17.9: Internal conversion background (col- 
ours) and signal (black dots) in the acoplanar mo- 
mentum - reconstructed mass plane for the phase 
I A detector configuration. 
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Figure 17.12: Tail of the internal conversion distri- 
bution overlaid with signal at different branching 
ratios for the phase IA detector. The resolution 
for the internal conversion decays was taken from 
30 000 simulated signal decays. 



Figure 17.10: Internal conversion background 
(colours) and signal (black dots) in the acoplanar 
momentum - reconstructed mass plane for the 
phase IB detector configuration. 
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Figure 17.11: Internal conversion background 
(colours) and signal (black dots) in the acoplanar 
momentum - reconstructed mass plane for the 
phase II detector configuration. 
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Figure 17.13: Tail of the internal conversion distri- 
bution overlaid with signal at different branching 
ratios for the phase IB detector. The resolution 
for the internal conversion decays was taken from 
30 000 simulated signal decays. 
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Figure 17.14: Tail of the internal conversion distri- 
bution overlaid with signal at different branching 
ratios for the phase II detector. The resolution 
for the internal conversion decays was taken from 
30 000 simulated signal decays. 

three Michel decays 3 x (/i — > evv) and from co- 
incidences of a radiative muon decay with in- 
ternal conversion \i -4 eee with a Michel decay 
[i — > evv. Both types of backgrounds are sig- 
nificantly reduced by applying vertex and time 
requirements. The combinatorial Michel back- 
ground depends quadratically on the vertex and 
time resolution, whereas the coincidences /i — > eee 
x /i — > evv coincidence rate scale linearly with the 
vertex and the time resolution. 

The vertex resolution in the transverse direction 
can be represented by the distance of closest ap- 
proach when the track is extrapolated to the muon 
decay position. This values is about 185 um in all 
phases of the experiment and given in Table |17.2| 
For a target size of 10 cm in length and 2 cm in 
diameter a reduction factor of 1 • 10~ 4 can be de- 
rived per coincidence. 

The estimation of the reduction factor from the 
timing cut, relevant for phases lb and II, is difficult 
as the design for the time of flight system has not 
been finalised yet. Another difficulty comes from 
the fact that the timing resolutions of the scintil- 
lating fibre detector the scintillating tiles detector 



are expected to be different and that particles de- 
pending on their flight direction will be measured 
in both detectors or by the scintillating fibre de- 
tector only if staying in the central region. In the 
latter case the particles will be measured many 
times in several turns, what also leads to an in- 
crease in precision. Preliminary studies indicate 

that timing resolutions of 200-300 ps in the scin- 
tillating fibre detector and < 100 ps in the Scin- 
tillating Tiles Detector can be achieved, see sec- 
tions [IT] and [12] For the sensitivity calculation 
it is assumed that all particles will be measured 
with a time resolution of better than 250 ps and 
that reduction factors of 5 • 10~ 3 per coincidence 
can be obtained at a signal efficiency of 90 %. 

17.4 Results 

The signal efficiencies and the expected back- 
ground reduction factors of the different selection 
cuts are summarized in Table |17.3| Single event 
sensitivities are calculated, which are defined as 
the branching ratio at which one background event 
is expected. The sensitivity of the experiment 
scales with one over the number of muons on tar- 
get as long as the single event sensitivity is not 
reached and with one over the square root of the 
number of muons on target after. It can be seen 
that the single event sensitivity at phase IB is 
higher than at phase II. This comes from the fact 
that the combinatorial background increases at 
higher muon rates. However, the time to reach 
the aimed sensitivity of 1 • 10~ 16 takes in phase IB 
ten years of running whereas at phase II this sens- 
itivity can be reached in 180 days because of the 
higher muon rate, see Figure |17.15| Although the 
expected sensitivity at phase IA, even without the 
time of flight system, is quite high with 4 • 10~ 16 
it would take a long time of running to reach this 
sensitivity. However, about one month of run- 
ning is sufficient to reach a sensitivity of 1 • 10~ 13 , 
which is factor ten smaller than the current exper- 
imental bound. 
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Phase IA 


Phase IB 


Phase II 


backgrounds: 








Michel 


o 


< 2.5 • 10~ 18 


5 • 10~ 18 


fi — > eeeuu 


1 • 10~ 16 


l'-io- 17 


1 • io- 17 


fi — > eeevv and accidental Michel 





< 2.5 • 10~ 21 


7.5 • 10~ 18 


Total Background 


1 • 10~ 16 


l-io- 17 


2.3 • 10~ 17 


Signal: 








m i i_ j ■ l i i_ ■ rr* 

irack reconstruction and selection efficiency 


26% 


39% 


38% 


Kinematic cut (2a) 


95 % 


95 % 


95 % 


Vertex efficiency (2.5er) 2 
Timing efficiency (2a) 2 


98% 


98% 


98% 




90% 


90% 


Total efficiency 


24% 


33% 


32% 


Sensitivity: 








Single event sensitivity 


4- 10~ 16 


3 • 10~ 17 


7- 10~ 17 


muons on target rate (Hz) 


2 • 10 7 


1 • 10 s 


2 • 10 9 


running days to reach 1 • 10 5 


2600 


350 


18 


running days to reach 1 • 10 -16 




3500 


180 


running days to reach single event sensitivity 


6500 


11700 


260 



Table 17.3: Signal efficiencies and estimated background reduction factors of the discussed selection cuts, 
which are used to determine the single event sensitivities for the three phases of the experiment. For the 
given muon on target rates running times are calculated to reach the given single event sensitivities. 
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Figure 17.15: Projected sensitivity and projected limit (90% CL) of the MuSe experiment for the three 
construction phases IA (red), IB (magenta) and II (blue) as function of the running time. The curves 
were obtained from the corresponding numbers in Table 17.3 For comparison also the (90% CL) obtained 
from the SINDRUM experiment is indicated (black dashed line). 
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The Mu3e Collaboration 
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The Institutes in Mu3e 



18.1 Responsibilities 

The responsibilities of the institutes participating 
in Mu3e are outlined in Table 118.11 This list con- 
cerns the phase I experiment and does not imply 
funding responsibilities. 

18.2 Collaborators 

18.2.1 University of Geneva 

Currently working on Mu3e: 

Alain Blondel Professor 

Alessandro Bravar Professor 

Martin Pohl Professor 

Antoaneta Damyanova Master student 

Planned positions: 

1 Postdoc 

1 Ph.D. candidate 

18.2.2 University of Zurich 

Currently working on Mu3e: 
Ueli Straumann Professor 
Peter Robmann Senior scientist 
Roman Gredig Ph.D. candidate 

18.2.3 Paul Scherrer Institut 

Currently working on Mu3e: 

Felix Berg Ph.D. candidate 

Malte Hildebrandt Senior scientist 

Peter-Raymond Kettle Senior scientist 

Angela Papa Senior scientist 



Stefan Ritt Senior scientist 

Alexey Stoykov Senior scientist 

Funded open positions: 

1 Postdoc 

1 Ph.D. candidate 

18.2.4 ETH Zurich 

Currently working on Mu3e: 
Christophorus Grab Professor 
Giinther Dissertori Professor 
Rainer Wallny Professor 
Planned positions: 
1 Ph.D. candidate 

18.2.5 University of Heidelberg 

Physikalisches Institut (PI) 

Currently working on Mu3e: 

Andre Schbning Professor 

Dirk Wiedner Senior scientist 

Sebastian Bachmann Senior scientist 

Niklaus Berger Junior research group leader 

Bernd Windelband Mechanical engineer 

Moritz Kiehn Ph.D. candidate 

Kevin Stumpf Mechanical engineer 

Raphael Philipp Master student 

Funded open positions: 

1 Postdoc 

4 Ph.D. candidates 

Planned positions: 

1 Postdoc 

1 Ph.D. candidate 
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Kirchhoff Institut fur Physik (KIP) 

Currently working on Mu3e: 

Hans-Christian Schultz-Coulon Professor 

Wei Shen Postdoc 

Patrick Eckert Ph.D. candidate 

Carlo Licciulli Master student 

Tobias Hartwig Bachelor student 

Planned positions: 

1 Postdoc 

1 Ph.D. candidate 

Zentralinstitut fur Technische Informatik 

Currently working on Mu3e: 
Peter Fischer Professor 
Ivan Peric Senior scientist 



Component 


Responsible 


Beam 


PSI 




Magnet 


UGS 






Heidelberg 


PI 


Target 


Heidelberg 


PI 




PSI 




Pixel chip 


ZITI Mannheim 




Heidelberg PI 




Heidelberg 


NB 


Pixel detector 


Heidelberg 


PI 


Mechanics and cooling 


Heidelberg 


PI 


Pixel on-detector electronics 


Heidelberg 


NB 




Heidelberg 


PI 


Read-out electronics 


Heidelberg 


NB 


i ibre detector 


Geneva 






Zurich 






ETHZ 






PSI 




Tile detector 


Heidelberg 


KIP 


Timing electronics 


PSI 






Heidelberg 


KIP 




Geneva 






Zurich 






ETHZ 




Filter farm 


Heidelberg 


NB 


Slow control 


PSI 




Infrastructure 


PSI 





Table 18.1: Institutional interests and responsib- 
ilities in phase I Mu3e. NB denotes the Emmy 
Noether junior reserch group at the Physikalisches 
Institut led by N. Berger. UGS is an external en- 
gineering consultancy tasked with magnet design 
and procurement. Note that especially for large 
common items such as the magnet, design and 
construction responisbility does not imply fund- 
ing responsibility 
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Schedule 



19.1 Phase I Schedule 

Our goal is to take first data in 2015, with a pixel- 
only detector. Correspondingly, the timing crit- 
ical items are the central pixel detectors and the 
magnet. Prototyping on the first is well underway 
and a first instrumented assembly corresponding 
to the vertex layers should be operational by late 
2013. Magnet design has also started, however the 
tendering process will consume about half a year; 
we hope the actual magnet construction will be 
under way early in 2014. 

For a second, longer run in 2016, we plan to in- 
tegrate the fibre detector as well as the inner recurl 
stations complete with timing tiles. As the outer 
recurl stations are identical to the inner ones, we 



plan to maintain the production capacities and 
manufacture them during 2016. 

Depending on the outcome of the 2016 run and 
the schedule for the high intensity beam line, 2017 
could either be spend for phase II preparations or 
with an extended run at the existing beamlinc. 

19.2 Phase II Schedule 

The phase II schedule is contingent on the sched- 
ule for the high intensity beamline. By 2017, we 
can have all required detector components ready. 
If the phase I experience shows that major parts of 
the detector will have to be rebuilt, earliest phase 
II startup would slip to 2018. 
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Cost Estimates 



Item 


Estimated cost 




(KCHF) 


Beam line 


< 50 


Target 


< 10 


Pixel detector (central) 


400 


Pixel detector (2 recurl stations) 


200 


Scintillating fibres 


200 


Tile detector (2 stations) 


110 


Mechanics and cooling 


200 


DAQ pixel detector 


50 


DAQ fibres 


200 


DAQ tiles 


190 


Central DAQ 


100 


Filter farm 


50 


Slow control 


50 


Infrastructure 


100 


Total 


1910 



Table 20.1: Estimated costs of the phase I experi- 
ment. The magnet (estimated at 500 KCHF) will 
be loaned to Mu3e by PI Heidelberg for the run 
time of the experiment. 



Item 


Estimated cost 




(KCHF) 


Beam line 


u.a. 


Target 


< 10 


Pixel detector (upgrade central) 


200 


Pixel detector (2 recurl stations) 


200 


Scintillating fibres 


200 


Tile detector (2 stations) 


110 


Mechanics and cooling 


100 


DAQ pixel detector 


50 


DAQ fibres 


200 


DAQ tiles 


190 


Central DAQ 


100 


Filter farm 


50 


Slow control 


50 


Infrastructure 


u.a. 


Total 


1460 



Table 20.2: Estimated costs of the phase II experi- 
ment. Items marked u.a. are under assesment; the 
new beam line required will be by far the most ex- 
pensive item but also benefit many other users. 
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Appendix 



A.l Mu3e theses 

Several theses involving the Mu3e project have 
been completed: 

• M. Kichn, diploma thesis Track Fitting with 
Broken Lines for the MuSe Experiment, 
Heidelberg University, 2012. 

• M. Zimmcrmann, bachelor thesis Cooling 
with Gaseous Helium for the MuSe Experi- 
ment, Heidelberg University, 2012. 

• H. Augustin, bachelor thesis Charakterisier- 
ung von HV-MAPS, Heidelberg University, 
2012. 

• A.-K. Perrevoort, master thesis Characterisa- 
tion of High Voltage Monolithic Active Pixel 
Sensors for the MuSe Experiment, Heidelberg 
University, 2012. 

These theses are available in port- 
able document format (PDF) from 
http : //www.psi . ch/mu3e/documents. Sev- 
eral more theses are ongoing (thus titles are 
preliminary) : 



• R. Gredig, doctoral thesis Fibre Tracker 
for the MuSe Experiment, Zurich University, 
started 2012. 

• M. Kiehn, doctoral thesis Track Reconstruc- 
tion for the MuSe Experiment, Heidelberg 
University, started 2012. 

• A. Damyanova, master thesis Fibre and SiPM 
characterization for the MuSe Experiment, 
Geneva University, started 2012. 

• R. Philipp, master thesis Tests of High 
Voltage Monolithic Active Pixel Sensors for 
the MuSe Experiment, Heidelberg University, 
started 2012. 
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